ARTIFICIAL  INHOMOGENEOUS 
TAPERED  IMPEDANCE  SHEET 
CHARACTERIZATION  AND  APPLICATIONS 


THESIS 


Jeffrey  P.  Massman,  Lieutenant,  USAF 
AFIT  /  GE/EN  G/ 10- 16 


DEPARTMENT  OF  THE  AIR  FORCE 
AIR  UNIVERSITY 

AIR  FORCE  INSTITUTE  OF  TECHNOLOGY 


Wright- Patterson  Air  Force  Base,  Ohio 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED. 


The  views  expressed  in  this  thesis  are  those  of  the  author  and  do  not  reflect  the 
official  policy  or  position  of  the  United  States  Air  Force,  Department  of  Defense,  or 
the  United  States  Government. 


AFIT/GE/ENG/10-16 


ARTIFICIAL  INHOMOGENEOUS  TAPERED  IMPEDANCE  SHEET 
CHARACTERIZATION  AND  APPLICATIONS 


THESIS 


Presented  to  the  Faculty 

Department  of  Electrical  and  Computer  Engineering 
Graduate  School  of  Engineering  and  Management 
Air  Force  Institute  of  Technology 
Air  University 

Air  Education  and  Training  Command 
in  Partial  Fulfillment  of  the  Requirements  for  the 
Degree  of  Master  of  Science  in  Electrical  Engineering 

Jeffrey  P.  Massman,  BSEE 
Lieutenant,  USAF 

March  2010 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED. 


AFIT/GE/ENG/10-16 


ARTIFICIAL  INHOMOGENEOUS  TAPERED  IMPEDANCE  SHEET 
CHARACTERIZATION  AND  APPLICATIONS 


Jeffrey  P.  Massman,  BSEE 
Lieutenant,  USAF 


Approved: 


.  Havrilla  (Chairman) 


rO  fj . 

Andrew  J.  Terzuoli  (Number) 


Cl  Moa.  A# i o 

Date 

7  MAR  Z01D 

Date 

9  M<yr  2fi/Q 


Date 


AFIT/GE/ENG/10-16 


Abstract 

The  ability  to  manipulate  a  complex  impedance  sheet  about  a  PEC  or  general 
dielectric  cylinder  of  arbitrary  diameter  to  control  the  overall  scattering  characteristics 
at  particular  angles  for  shielding  or  scattering  applications  is  presented.  It  is  shown 
that  with  the  capability  to  readily  control  both  the  real  and  imaginary  components 
of  a  percolated  lossy  film’s  effective  impedance,  specific  impedance  distributions  can 
be  obtained  to  effectively  direct  the  overall  nulls  in  the  bistatic  scattering  pattern. 
The  exact  approach  only  requires  negligible  truncation  error  of  the  matrix  solution 
and  can  be  applied  to  cylinders  of  arbitrary  size.  This  approach  relies  on  general 
impedance  sheet  boundary  conditions  in  addition  to  a  Fourier  series  representation  of 
the  impedance  sheet  distribution,  with  the  utilization  of  Wronskian  and  recurrence 
relations,  to  produce  a  computationally  efficient  matrix  solution  for  the  scattered  field 
expansion  coefficient  an.  This  scattering  coefficient  an  can  then  be  included  in  the 
expression  for  a  two-dimensional  scattering  width,  allowing  for  the  computation  of 
scattering  patterns  of  arbitrary  diameter  cylinder  impedance  sheet  distributions.  The 
ability  to  adapt  sheet  impedances  to  effectively  control  the  overall  nulls  in  the  bistatic 
scattering  pattern  will  be  demonstrated,  with  practically-implemented  impedance 
characteristics.  Lastly,  a  novel  rectangular  X-band  waveguide  system  will  be  utilized 
to  characterize  the  electromagnetic  properties  of  the  ^-dependent  impedance  sheet. 
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ARTIFICIAL  INHOMOGENEOUS  TAPERED  IMPEDANCE  SHEET 
CHARACTERIZATION  AND  APPLICATIONS 

I.  Introduction 

In  the  past  decade  electromagnetic  compatibility  (EMC)  issues  have  increased  in 
stature  as  low  emission  and  high  immunity  to  interference  have  emerged  as  essen¬ 
tial  criteria  for  characterizing  overall  integrated  circuit  (IC)  performance.  Achieving 
low  emissions  and  increased  interference  immunity,  however,  can  by  problematic  as 
technological  trends  have  a  propensity  to  act  in  an  opposing  manner.  For  instance, 
advances  in  process  integration,  higher  switching  speeds,  and  more  complex  circuits 
generally  result  in  ICs  generating  an  increased  amount  of  parasitic  emissions.  Simi¬ 
larly,  immunity  to  electromagnetic  interference  (EMI)  tends  to  decrease  with  reduced 
supply  voltages  and  increasing  numbers  of  interfaces.  As  a  consequence,  the  EM- 
C/EMI  held  at  the  component  level  has  gained  much  attention  to  abet  with  IC  design 
in  order  to  achieve  these  two  criterionfl], 

1.1  Problem  Statement 

Resistive  and  impedance  sheets  find  numerous  uses  in  EMC  and  EMI  applications. 
With  recent  advances  in  the  production  of  carbon  loaded  dielectric  films,  such  as  the 
ability  to  simultaneously  alter  the  real  and  imaginary  components  of  an  impedance 
sheet  [2],  novel  impedance  distributions  can  be  realized.  These  new  impedance  dis¬ 
tributions  can  allow  for  precision  scattering  and  shielding  control  and  manipulation, 
as  will  be  demonstrated  in  this  research.  Additionally,  recent  interest  in  the  scatter¬ 
ing  measurements  of  uniform  impedance  sheets  wrapped  about  a  polystyrene  foam 
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cylinder  was  reported  in  [3].  For  the  purposes  of  the  current  analysis,  the  problem 
is  intensified  by  the  impedance  sheet  distribution  having  a  general  0  dependency. 
While  there  exists  an  abundance  of  literature  investigating  resistive  and  impedance 
sheet  surfaces  and  scattering  [4]- [13],  this  study  introduces  an  analytic  approach  for 
calculating  the  scattering  characteristics  from  cylindrical  sheet  impedances  having  a 
0  dependency. 

Furthermore,  the  0  dependency  presents  an  additional  degree  of  complexity  for 
electromagnetic  material  characterization  purposes.  The  process  of  determining  the 
complex  permittivity  and  permeability  of  a  material  typically  requires  the  sample 
to  be  linear,  homogeneous,  and  isotropic.  The  introduction  of  the  0  dependency 
explicitly  violates  the  homogeneous  requirement.  To  mitigate  inhomogeneous  effects 
a  reduced  aperture  waveguide  system  will  be  presented.  This  stepped  waveguide 
technique  utilizes  a  modal  analysis  to  account  for  the  excited  higher-order  modes, 
effectively  allowing  smaller  cross-sections  of  the  0-varying  impedance  sheet  to  be 
characterized  [14;  15]. 


Figure  1.  Experimental  setup  variations  of  the  full  X-band  plates,  stepped  aperture 
plates  of  0.1,  0.2,  0.3  inches,  and  stepped  aperture  sheets  of  0.1,  0.2,  0.3  inches. 
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1.2  Limitations 


Despite  employing  an  analytical  approach  to  determining  the  scattering  and  mate¬ 
rial  characteristics  from  the  0- varying  impedance  sheet,  several  practical  aspects  must 
be  understood.  For  instance,  the  novel  impedance  sheet  distributions  necessary  to  at¬ 
tain  precision  scattering  and  shielding  control  and  manipulation  require  ideal  material 
values.  Introducing  realistic  values  diminishes  the  distribution  effectiveness,  such  as 
frequency  dependence.  Additionally,  the  reduced  aperture  waveguide  material  char¬ 
acterization  technique  limits  the  amount  of  energy  transmitted  into  the  impedance 
sheet  sample,  thereby  increasing  the  material  characterization  uncertainty.  To  exac¬ 
erbate  the  problem  an  impedance  sheet  distribution  containing  high-loss  aspects,  to 
include  low  frequency  measurements,  can  lead  to  total  reflection  with  transmission 
measurements  near  the  noise  floor. 

1.3  Scope 

This  research  will  include  a  theoretical  development  of  the  TM2  and  TE2  held  sets 
for  a  cylinder  comprised  of  a  (^-dependent  impedance  sheet.  General  impedance  sheet 
boundary  conditions  will  be  applied  in  addition  to  a  Fourier  series  representation 
of  the  impedance  sheet  distribution.  The  utilization  of  Wronskian  and  recurrence 
relations  will  help  to  produce  a  computationally  efficient  matrix  solution  for  the  scat¬ 
tered  held  expansion  coefficient  an.  This  scattering  coefficient  an  will  then  be  included 
in  the  development  of  a  two-dimensional  scattering  width,  resulting  in  the  scatter¬ 
ing  patterns  of  arbitrary  diameter  cylinder  impedance  sheet  distributions.  Then, 
the  ability  to  adapt  sheet  impedances  to  effectively  control  the  overall  nulls  in  the 
bistatic  scattering  pattern  will  be  demonstrated.  Lastly,  the  electromagnetic  mate¬ 
rial  characterization  aspect  will  focus  on  X-band  frequencies  (8.2  to  12.4  GHz)  within 
a  rectangular  waveguide  system,  where  a  Newton-Raphson  root  search  algorithm 


3 


will  iteratively  compute  the  complex  permittivity  of  the  non-magnetic  ^-dependent 
impedance  sheet. 

1.4  Thesis  Organization 

Chapter  2  provides  the  necessary  theoretical  background  beginning  with  constitu¬ 
tive  parameters  and  relations  and  Maxwell’s  equations  in  time-harmonic  form,  deriv¬ 
ing  the  impedance  boundary  condition,  and  leading  to  the  computationally-efhcient 
incident  TM2  and  TE2  held  sets.  Chapter  3  derives  an  analytic  approach  for  cal¬ 
culating  the  scattering  characteristics  from  cylindrical  sheet  impedances  having  a  0 
dependency  with  free-space  filled,  partially- filled  perfect  electric  conductor  (PEC), 
and  general  dielectric  filled  cases.  Chapter  4  illustrates  how  impedance  distributions 
can  allow  for  precision  scattering  and  shielding  control  and  manipulation  with  exper¬ 
imentally  measured  complex  material  parameters.  Finally,  Chapter  5  provides  a  brief 
set  of  conclusions  and  recommendations  for  future  research. 
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II.  Background  Electromagnetic  Theory 


The  study  of  electromagnetic  field  theory  is  concerned  with  charges  in  a  general 
state  of  motion  that  produce  electric  and  magnetic  fields.  The  mathematical  model 
that  characterizes  the  relationship  between  the  electromagnetic  fields  and  the  sources 
which  maintain  these  fields  is  known  as  Maxwell’s  equations.  For  this  research,  all 
electromagnetic  phenomena  will  be  viewed  from  the  macroscopic  perspective,  where 
the  space-time  average  of  sources  and  fields  over  regions  and  periods  are  large  in 
comparison  to  atomic  dimensions.  Additionally,  the  rationalized  rnksc  system  of 
units  is  employed.  Finally,  all  theoretical  development  revealed  in  this  chapter  can 
be  readily  followed  via  [16]- [21], 

2.1  Constitutive  Parameters  and  Relations 

Materials  contain  charged  particles,  such  as  electrons  and  protons,  that  necessarily 
interact  with  externally  applied  electric  and  magnetic  fields  based  on  their  respective 
atomic  structures.  This  interaction  is  physically  characterized  through  constitutive 
parameters  which  define  and  mathematically  model  the  relationship  between  charges 
and  fields.  These  constitutive  parameters  are  known  as  permittivity,  permeability,  and 
conductivity.  The  extent  of  this  research  assumes  the  complex  cylindrical  impedance 
sheet  is  immersed  in  simple  medium  (linear,  homogeneous,  and  isotropic),  as  com¬ 
monly  defined  as 


Mc  =  0, 

(1) 

Jc  (7  E , 

(2) 

D  =  eE, 

(3) 

B  =  /iH, 

(4) 
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where  Mc  is  the  magnetic  conduction  current  density  (volts  /  square  meter),  Jc  is 
the  electric  conduction  current  density  (Amperes  /  square  meter),  er  is  the  electric 
conductivity,  D  is  the  electric  flux  density  (Coulombs  /  square  meter),  E  is  the  electric 
field  intensity  (Volts  /  meter),  e  is  the  complex  permittivity  (Farads  /  meter),  B  is 
the  magnetic  flux  density  (Webers  /  square  meter),  H  is  the  magnetic  field  intensity 
(Amperes  /  meter),  and  /i  is  the  complex  permeability  (Henries  /  meter). 

For  simple  medium,  the  permittivity  parameter  e  defines  the  susceptibility  of  a 
material  to  being  polarized  by  an  external  electric  field.  Likewise,  the  permeability 
parameter  //  defines  the  susceptibility  of  a  material  to  being  magnetized  by  an  external 
magnetic  field.  The  conductivity  parameter  establishes  the  susceptibility  of  a  material 
to  sustain  an  electric  (a)  or  fictitious  magnetic  ( am )  conduction  current  given  an 
applied  electric  or  magnetic  field.  Linearity  relates  the  induced  sources  to  the  field 
strength;  homogeneity  describes  the  spatial  variation;  and  isotropic  means  the  induced 
currents  are  aligned  with  the  applied  fields. 

2.2  Time-Harmonic  Electromagnetic  Fields 

The  theory  and  analysis  presented  in  this  research  focuses  solely  on  electromag¬ 
netic  phenomena  that  vary  sinusoidally  in  time.  This  steady-state  sinusoidal  time 
dependence  leads  to  a  mathematically  simplified  phasor  notation  of  Maxwell’s  equa¬ 
tions,  such  that 


V  x  E(r,u>)  =  —M(f,  to)  —  juB(r,u),  (5) 

V  x  H(f,u)  =  J(r,  u)  +  juD(r,u),  (6) 

V  ■  D{f ,  u>)  =  qev(r  ,cj),  (7) 

V  •  B(f,  u)  =  qmv(r,  u),  (8) 
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where  qev  is  the  electric  charge  density,  qmv  is  the  magnetic  charge  density,  f  is  the 
position  vector,  and  uj  is  the  angular  excitation  frequency.  Substituting  (l)-(4)  into 
(5)-(8)  leads  to  the  following  (if  source  free) 


V  x  E(r,ui)  =  —jufiH(r,uj), 

(9) 

V  x  H(r,u>)  =  oE  +  jcjeE(r,Lu), 

(10) 

ri  Br*  \  Qev(r,uj) 

V  •  E{r ,  lo)  — 

(H) 

Y7  \  9mv(r,u;) 

V-H(r,u)  = 

(12) 

h 

which  represents  the  time-harmonic  Maxwell’s  equations  in  simple  media. 

2.2.1  Boundary  Conditions. 

The  macroscopic  electromagnetic  held  behavior  across  a  boundary  can  be  discon¬ 
tinuous.  The  boundary  conditions  for  a  generalized  impedance  sheet  will  be  derived 
from  Maxwell’s  equations.  Boundary  value  problems  explored  in  later  sections  will 
be  solved  via  imposing  the  resulting  conditions.  A  more  complete  understanding  of 
the  boundary  conditions  can  illuminate  the  limitations  and  overall  usefulness  of  the 
implemented  technique. 


where  E\  Hl  are  the  incident  electric  and  magnetic  fields,  ES,HS  are  the  scattered 
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electric  and  magnetic  fields,  Ef,  H*  are  the  transmitted  electric  and  magnetic  fields, 
eo  is  the  free  space  permittivity,  po  is  the  free  space  permeability,  and  A  is  the  elec¬ 
tromagnetic  field’s  wavelength.  Adapting  (9)  and  (10)  to  the  boundary,  as  shown 
in  Figure  2,  leads  to  the  following  formulation  (assuming  the  impedance  sheet  is 
non- magnetic) 


V  x  E  =  —juj/doH, 
V  x  H  =  <t E  +  jcoeE. 


(13) 

(14) 


Next,  integrating  (13)  throughout  the  enclosed  volume  shown  in  Figure  2  and  allowing 
the  volume  to  approach  zero  produces 


lim  /  V  x  EdV  =  —juJUo  dm  /  HdV. 
AV^O  J Av  AV^O  JAV 


(15) 


If  H  is  bounded  as  AV 


0,  then  the  right-hand  side  of  (15)  goes  to  zero,  mathe¬ 


matically  shown  as 


— juj/j,0  ^linio  j  HdV  =  —jun0H  AV  =  0.  (16) 

Then,  the  Curl  Theorem  may  be  used  to  convert  the  volume  integral  of  the  remaining 
portion  from  (15)  to  a  surface  integral,  such  that 


n  x  EdS  =  0. 


AS 


(17) 


If  the  radius  R  is  much  larger  than  the  electromagnetic  wavelength,  then  the  expres- 


sion  can  be  approximated  such  that 


n  x  EdS  ~  [n  x  E+  —  n  x  E  +  ns  x  E)AS  =  0, 


AS 


(18) 


where  E±  =  E{R±,^>)  and  n  is  the  unit  outward  normal.  Here,  large  contrast  is  the 
key,  depicted  as 


Figure  3.  High  contrast  electric  field  at  the  impedance  boundary. 


where  E  is  the  tangential  electric  field.  The  large  contrast  results  in  hs  x  E  —  0  and 
leads  to  the  following  boundary  condition 


n  x  E+ 


h  x  E 


(19) 


The  second  boundary  condition  may  be  found  by  first  adding  a  zero  sum  pair  of 
terms  to  the  right-hand  side  of  (14)  resulting  in 


aE  +  jueE , 

(20) 

crE  +  jueE  —  jue^E  +  jueoE, 

(21) 

aeE  +  jue0E, 

(22) 

where  ae  =  a  +  ju(e  —  eo).  This  leads  to  the  notion  of  effective  or  equivalent  conduc¬ 
tivity  immersed  in  free-space.  Then,  in  an  analogous  process,  integrating  (22)  as  the 
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volume  approaches  zero  results  in 


lim  [  V  x  HdV  =  [  JedV  +  jue0  [  EdV,  (23) 

A^o  JAV  JAV  JAV 

where  the  bounded  equivalent  conduction  current  is  defined  as  Je  =  aeE.  With  E 
bounded  the  second  right-hand  term  of  (23)  must  go  to  zero  as  given  by 

jue o  f  EdV  =  0.  (24) 

JAV 

Again,  the  Curl  Theorem  may  be  applied  to  the  effect  that 

n  x  (. H+  -  H~)AS  =  Je AV,  (25) 

where  the  volume  differential  may  be  expanded  in  terms  of  a  surface  differential, 
JeAV  =  JedAS  =  JesAS.  This  results  in  (25)  reducing  such  that 

n  x  (H+-H~)  =  Jes.  (26) 

If  it  is  then  assumed  that  the  resistive-air  interfacial  contrast  is  sufficient,  then  an 
equivalent  impedance  can  be  introduced  as 


Jes  Jed 

(27) 

(&edJtang)d 

(28) 

1 

(29) 

Etang) 

^ e 

where  Etang  is  the  tangential  electric  field  and  the  equivalent  impedance  is  defined 
as  Ze  =  Aj.  Substituting  (29)  into  (26)  leads  to  the  following  boundary  condition 
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formulation 


n  x  (n  x  E^)  =  —  Zen  x  ( H+  —  H  ).  (30) 

As  a  set  of  equations  both  (19)  and  (30)  constitute  generalized  approximate  bound¬ 
ary  conditions  for  a  thin  non-magnetic  impedance  structure.  To  summarize,  the  gen¬ 
eralized  impedance  boundary  conditions  are  presented  as 

n  x  E+  —  h  x  E~  (31) 

n  x  {n  x  E±)  =  —Zeh  x  ( H+  —  H~ ).  (32) 

In  the  physical  sense,  (31)  implies  that  the  electric  held  is  continuous  across  the 
impedance  sheet,  while  (32)  indicates  that  the  magnetic  held  is  discontinuous  across 
the  impedance  boundary  by  the  amount  of  electric  current  induced  in  the  impedance 
sheet.  This  can  be  readily  perceived  with  an  equivalent  circuit  theory  model. 


Figure  4.  Impedance  sheet  boundary  condition  circuit  analogy. 

Figure  4  models  the  electrically  thin  impedance  sheet  with  a  lumped  impedance 
element  (Ze)  in  a  parallel  configuration.  The  electric  holds  on  either  side  of  the 
boundary  are  represented  by  the  voltages  V\  and  V2  and  are  equal  in  the  parallel 
conhguration.  The  magnetic  helds,  represented  by  the  currents,  illustrate  how  the 


11 


current  from  V\  to  V2  is  discontinuous  by  the  amount  of  shunt  current  that  flows 
through  the  parallel  impedance  element  Ze. 

2.3  Scattering  by  Cylinders 

The  extent  of  this  research  will  focus  on  modal  analysis  for  held  solutions.  In 
particular  the  scattering  characteristics  from  cylindrical  sheet  impedances  having  a  (j) 
dependency  will  be  explored.  The  development  utilizes  a  two-dimensional  approach 
in  order  to  obtain  a  meaningful  solution  with  an  accompanied  physical  intuition  not 
readily  apparent  or  attainable  from  a  three-dimensional  development.  For  practical 
implementation  an  approximate  2D-to-3D  expansion  may  be  employed 

2  l2 

&3D  ~  C2D— (33) 

where  L  is  the  physical  height  of  the  cylinder,  which  is  assumed  to  be  suitably  located 
for  broadside  illumination. 

2.3.1  Electromagnetic  Scattering  Methodologies. 

Electromagnetic  scattering  concerns  the  phenomena  of  an  impressed  source  exci¬ 
tation  inducing  charges  and  currents  on  a  structure,  which  then  re-emits  the  energy 
in  the  form  of  electromagnetic  waves.  The  general  approach  to  a  solution  is  to  deter¬ 
mine  the  scattered  fields  Es,  Hs  since  they  can  then  be  utilized  for  RCS,  radiation 
pattern,  and  impedance  calculations.  With  many  existing  methods  for  determining 
the  scattered  fields,  the  particular  choice  method  typically  depends  upon  the  geomet¬ 
rical  shape  and  electrical  size  of  the  scattering  body.  These  techniques  are  classified 
analytic,  low-frequency,  or  high-frequency  methods.  For  this  research,  an  analytic 
method  based  upon  a  series  expansion  will  be  implemented.  Essentially,  the  modal 
expansion  technique  is  a  field-based  approach  such  that  the  fields  are  expanded  in 
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infinite  series  of  modes  where  each  mode  is  a  solution  to  the  wave  equation.  Bound¬ 
ary  condition  enforcement  results  in  infinite-series  expansion  coefficients  which  can 
then  be  applied  to  achieve  an  exact  infinite-series  solution  for  the  scattered  fields. 
This  technique  has  inherent  advantages  and  disadvantages.  It  is  exact  and  rigorous 
but  may  be  both  computationally  intensive  and  slow  to  converge  if  the  scattering 
structure  is  much  larger  than  the  electromagnetic  excitation  frequency. 

2.3.2  Wave  Decoupling  due  to  z-invariance. 

With  a  two-dimensional  approach  the  concept  of  invariance  may  be  introduced. 
This  requires  two  additional  assumptions  placed  upon  the  problem  setup.  First, 
it  necessitates  that  an  excited  electromagnetic  wave  be  normally  incident  upon  the 
cylinder.  For  the  development  of  this  analysis,  the  plane  wave  will  be  expressed  as 

E\p)  =  zEl(x)  =  zE'0ze-*kox>  (34) 

H\p)  =  zHt(x)  =  zHl0ze~jkoX.  (35) 

Second,  the  cylinder  will  be  assumed  to  have  a  length  infinite  in  extend  in  the  ,2- 
direction,  as  shown  in  Figure  5. 


Figure  5.  2D  z-invariance  physical  structure. 
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With  both  the  geometrical  structure  and  the  excitation  being  z- invariant,  the  incident 
and  scattered  waves  will  decouple  into  two  independent  field  sets,  TEZ  and  TM~. 
This  is  mathematically  demonstrated  by  starting  with  Maxwell’s  equations  and  the 
continuity  equations,  such  that 


V  x  E  =  —M  —  jujfioH , 

(36) 

V  x  H  —  J  - b  jco6()E , 

(37) 

V  *  -D  Qevi 

(38) 

V  *  -B  Qmvi 

(39) 

Due  to  the  inherent  ^-invariance  of  the  problem,  the  5  component  of  the  differential 
operator  is  eliminated,  fz  =  0.  This  is  not  to  say  that  a  ^-polarized  electric  field 
cannot  exist,  but  rather  that  if  it  does  then  Ez  must  be  constant  for  all  This 
causes  the  differential  operator  (for  cylindrical  coordinates)  to  reduce  such  that 


(40) 


where  the  t  subscript  has  been  adopted  to  indicate  a  transverse  operator.  Applying 
(40)  to  (36)-(39)  and  the  continuity  equations  results  in 


Vt  x  (Et  +  zEz )  =  — ( Mt  +  zMz)  —  +  zHz ),  (41) 

Vi  x  (Ht  +  zHz )  =  ( Jt  +  zJz )  +  jue0{Et  +  zEz ),  (42) 

Vi  •  (Et  +  zEz)  =  — ,  Vi  •  (Jt  +  zJz )  =  —jujqeVi  (43) 

Co 

Vi  ■  (Ht  +  zHz )  =  ^,  Vi  •  (Mt  +  zMz )  =  -juqmv.  (44) 

ho 
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Expanding  and  grouping  terms  in  each  equation  leads  to 


V(X£,  +  V,x  zEz  =  -Mt  -  juJUoHt  ~  zMz  -  zjuii0Hz,  (45) 

V*  x  Ht  +  Vi  x  zHz  =  Jt  +  jue0Et  +  zJz  +  zjuj6oEz,  (46) 

VfEt  =  —,  VfJt  =  -juqev,  (47) 

eo 

Vt-Ht  =  •  Mt  =  -juqmv.  (48) 

ho 


Next  equating  transverse  and  longitudinal  terms  in  the  curl  equations  and  identifying 
transverse  component  terms  in  the  divergence  and  continuity  equations  leads  to  the 
following  decoupled,  independent  sets  of  equations 


TE~  =  < 


TM~  =  < 


VtxEt  =  —zMz  -  zjun0Hz 
Vi  x  zHz  =  Jt  +  jue0Et  , 
y.  ■  Et  =  ^ 

1  L  eo  jueo 

Vi  x  Ht  =  zJz  +  zju;e0Ez 
Vi  x  zEz  =  -Mt  -  junoHt  > 

.  jj  _  jmii  __  -Vt-Mt 
t  t  Mo  jtoMO 


(49) 


(50) 


where  the  Duality  theorem  principle  holds  and  can  be  applied  to  interchange  between 
the  two  held  sets.  The  process  of  equating  the  terms  above  relied  on  the  inherent 
linear  independence  of  the  held  components. 


2.4  Cylindrical  Wave  Transformations  and  Theorems 

In  order  to  simplify  boundary  condition  enforcement  for  the  cylindrical  scatterers, 
various  cylindrical  wave  transformations  and  theorems  will  be  adopted.  Using  the 
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geometry  presented  in  Figure  5,  an  incident  plane  wave  can  be  represented  as 


E\p)  =  zEt(x)  =  zEl0ze^kopcos^,  (51) 

H\p)  =  zHt(x)  =  zHl0ze~jkopcoS(t>,  (52) 


if  E\  H‘  is  propagating  along  the  +x  direction,  for  the  TMZ  and  TEZ  helds  respect¬ 
fully.  For  the  cylindrical  wave  transformation,  representing  the  e_j7cpcos^  term  as 
Bessel  functions  will  take  into  account  the  radial  dependence  while  periodic  functions 
account  for  the  0-variation.  To  maintain  held  generality  the  origin  must  be  included, 
which  disallows  the  use  of  Bessel  functions  of  the  second  kind.  Additionally,  the 
periodicity  requirement  can  be  satisfied  via  a  term.  This  leads  to  the  following 
mathematical  representation 


OO 

e-jk0pcos<t>  =  anJn(k0p)e 

n=— oo 


(53) 


The  unknown  expansion  coefficient  an  can  then  be  found  via  integrating  both  sides 
of  (53)  by  the  operator  resulting  in 


r»  2n 


r*2n 


e~  jkoP  cos  <j>e~j 


m4d<j)=  ^  anJn(k0p)  /  e 


j(n—m. 


(54) 


where 


27 r  . .  .n  —  m 
0  . .  .n  7^  m 


(55) 


Ignoring  the  trivial  case  leads  to  the  following  formulation  (for  n  =  m) 


am2nJm(kp). 


(56) 
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Reindexing  (56)  and  applying  the  following  definition  and  identity 


j-n  r  27T 

J  „(kop)  =  —  /  e~jkopcos(pejn*d(j), 

2vr  Jo 

(57) 

Jn(^dP)  =  J— ro(  kop) 

(58) 

leads  to 


27 rj  nJn(k0p)  =  an27rJn(fc0p). 


(59) 


Thus,  an  =  j  n  and  the  cylindrical  wave  function  representation  of  the  plane  wave 
can  be  written  as 


e~jk0pcos<t>  =  ^  j-nJn(k0p)ejn*.  (60) 


It  then  follows  that  the  incident  electromagnetic  fields  for  the  TM2  and  TE2  case  can 
be  written  as 


OO 

E'(p)  =  zE‘0le~i =  zEl  V  (61) 

n=— oo 

OO 

Wifi)  =  zH'0ze~^  =  zW0,  V  (62) 

n=— oo 


which  can  be  manipulated  into  a  more  computationally  efficient  form  through  the  use 
of  certain  Bessel  function  properties,  as  shown  in  Appendix  A,  represented  as 


E\p)  =  zE'0ze-*kox  =  zEqz  enj~nJn(k0p)  cos  n0,  (63) 

n= 0 
oo 

H\p)  =  zW0ze~jkoX  =  zHl0zJ2tnrnUkop)cosn<f>,  (64) 

n= 0 
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where 


{1  . . .  n  =  0 

2  ...n^O 


(65) 


2.4.1  TM2  Polarization  Development. 

For  TM2  polarization,  the  incident  electric  held  expressed  in  (63)  can  be  used  to 
determine  the  incident  magnetic  held  components.  With  the  incident  electric  held  in  a 
convenient  cylindrical  wave-transformation  form,  the  previously  developed  z-invariant 
Faraday’s  law  may  be  applied  to  find  the  incident  magnetic  held.  Additionally  it  can 
be  assumed  that  the  TM2  incident  held  is  maintained  solely  by  a  physically  realizable 
impressed  ^-directed  current  density  Jl  with  Ml  =  0.  Thus,  the  incident  magnetic 
held  can  be  found  from  (50)  via 


-Vt  X  (zEi) 
jup  0 

(66) 

5  x  Vf Elz 

jun  o 

(67) 

1  .  ,JEi  2l5Ei. 

zx(p  z+(j)  2 

jkoVo  dp  p  d(p 

(68) 

1  .  SE*  21 5EI. 

zx(pk0  +0 

jk0r ?o  o(k0p)  p  d (p 

(69) 

.  1  dEi  !  ^  1  5Et 

Pjk0VoP  d(j)  9 jr]0d(k0p)' 

(70) 

The  magnetic  held  components  in  computationally  efficient  form  may  may  then  be 
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found  by  substituting  (63)  into  (TO),  such  that 


&  1 
p  jkoVop  <50 

(71) 

1  r  00 

=  -U  r  ,  [E0z  Y]  enj  nJn(W)cOS770] 

jkoVoP  0(p  ^ 

(72) 

npi  °° 

=  0  V  nj  nJn{kQp)  sin 770, 

JkoVop  “ 

(73) 

Hi  1  *El 

0  jVoS(k0p) 

(74) 

\  s  °° 

=  •  Xfl  ^  iE0z  Y]  enj  nJn(/h)P)cOST70] 

JVo  °{ko  P)  ^ 

(75) 

i  °° 

=  —  y,  enj~nJ(0A;op)  cos  7i0. 

(76) 

To  summarize,  the  TM2  incident  electromagnetic  field  components  expressed  in  com¬ 
putationally  efficient  form  are  shown  below 


OO 

El(p,  <P)  =  E0z  eni~njn(k)p)  COS  710, 
n= 0 

Hp(p ,  0)  =  2E°Z  jr  nj~n3n(k0p )  sin  7i0, 
jkoVoP  “ 

°° 

ms  no. 


(77) 

(78) 

(79) 


2.4.2  TES  Polarization  Development. 

For  TE2  polarization,  the  incident  magnetic  field  expressed  in  (64)  can  be  nsed  to 
determine  the  incident  electric  field  components.  With  the  incident  magnetic  field  in  a 
convenient  cylindrical  wave-transformation  form,  the  previously  developed  z-invariant 
Ampere’s  law  may  be  applied  to  find  the  incident  magnetic  field.  Additionally  it  can 
be  assumed  that  the  TE2  incident  field  is  maintained  solely  by  an  impressed  i-clirected 
current  density  M\  with  ,/)  =  0.  Thus,  the  incident  electric  field  can  be  found  from 
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(49)  via 


v,  x  (zHi) 

jU60 

(80) 

z  x  VtHlz 
jue  o 

(81) 

Vo  ,  ,A6Hi  26H* 

-^zx  (p-— -  +  ^, f 
jk0  p  dp  d(p 

(82) 

Vo  ,  t,k  5 Hi  2SHi. 

■i  zx(p  xd  ~\+0  xD 
jko  pd{k0p)  0(() 

(83) 

,  Vo  SHlg  7  Vo  SHI 

Pjk0p  d<p  j  S(kopy 

(84) 

The  magnetic  field  components  in  computationally  efficient  form  may  then  be  found 
by  substituting  (64)  into  (84),  such  that 


Vo  SHI 
jk0p  8(p 

(85) 

■l  AX  ^  "Jn(M  C0Sn$ 

jkopH  ^ 

(86) 

0z  V  nj  nJn(kop )  sin  n(f), 

(87) 

oo  w; 

j  d(kop) 

(88) 

■  r, ,  ,  [AS,  eni  nJn(^0d)  COS  U0] 

J  <KM  ^ 

(89) 

oo 

jHlozVo  ^  enj~n3'n(k0p )  cosn<£. 

(90) 

n= 0 

To  summarize,  the  TEZ  incident  electromagnetic  field  components  expressed  in  com- 
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putationally  efficient  form  are  shown  below 


oo 

H'z(p ,  0)  =  E0z  enj~nJn(koP)  COS  710, 

71— 0 

(91) 

Ep(pA)=  .,°  j^nj  "Jn(W)smn0, 

(92) 

OO 

El(p,  <t>)  =  jHLvo  e^'~njn(fcoP)  cos n(j). 

n= 0 

(93) 
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III.  Analytical  Methodology 


This  research  will  demonstrate  the  ability  to  adapt  sheet  impedances  to  effec¬ 
tively  control  the  overall  nulls  in  the  bistatic  scattering  pattern.  Utilizing  a  cylinder 
comprised  of  a  ^-dependent  impedance  sheet  provides  for  an  additional  degree  of 
complexity  beyond  uniform  impedance  sheets.  This  additional  complexity,  however, 
prevents  a  closed-from  solution  from  being  easily  obtained.  The  general  approach  to 
achieving  a  modal  solution  is  to  define  the  electromagnetic  fields  for  all  regions  in  a 
suitable  representation,  enforce  appropriate  boundary  conditions,  and  solve  for  the 
scattered  field  expansion  coefficient.  With  this  scheme  an  orthogonal  testing  operator 
is  typically  utilized,  after  enforcing  the  boundary  conditions,  to  identify  only  those 
electromagnetic  modes,  or  held  configurations,  which  satisfy  both  Maxwell’s  equa¬ 
tions  and  the  boundary  conditions.  This  leads  to  a  unique  and  physically  meaningful 
solution.  The  complexity  of  the  ^-dependent  impedance  sheet,  however,  precludes 
the  testing  operator  from  sifting  through  the  infinite  number  of  modes.  Hence,  a 
computationally  efficient  approach  will  be  utilized  and  demonstrated  in  the  following 
section. 

3.1  Sheet  Impedance  Distribution  Formulation 

Expanding  the  impedance  parameter  Zs((j>)  into  a  Fourier  series  will  allow  for  a 
computationally-efficient  formulation  to  identify  the  expansion  coefficients.  A  suitable 
and  physically  realizable  representation  for  the  impedance  distribution  about  the 
cylinder  to  achieve  this  outcome  is 

Q 

Zs{4> )  =  ^aq  cos (g0  +  (5q).  (94) 

q= o 

An  appropriate  testing  operator  to  later  assist  in  readily  determining  the  expansion 
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coefficients  after  boundary  condition  enforcement  is 


7 T 


cos  m<j){}d(j). 


(95) 


Applying  (95)  to  (94),  in  addition  to  a  cos  n(f>  term  extracted  in  advance  from  the 
computationally  efficient  TMZ  and  TE2  field  sets,  results  in 


Cmn  —  ~  /  Zs{(j))  cos m(f> cos  n(f>d(f>, 

7T 


(96) 


where  the  Cmn  label  has  been  adopted  to  signify  the  two-element  expansion  in  matrix 
form.  This  term  can  be  further  expanded  and  simplified  such  that 


Cmn  =  ~  /  Zs((j>)  cos m(j) cos n(j)d(j) 


7T 
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Q 


aq  cos {qcj)  +  (3q)  cos  rrup  cos  ruj)d(j) 


'  q= 0 

Q 


=  -J2ag  L  cos (q4>  +  j3q)  cos  mcj)  cos  n(j)d(j). 


q= o 


(97) 

(98) 

(99) 


Within  (99)  the  (/-indexed  term  can  be  re-written  utilizing  a  simple  cosine  identity  as 


cos  (q(f>  +  (3q)  =  cos  qcj)  cos  (3q  —  sin  qcf)  sin  (3q, 


(100) 


which  can  be  substituted  back  in  and  arranged  such  that 


Q 


Cmn  E  aq  /  [cos  qcj)  cos  f3q  cos  rri<p  cos  nti  —  sin  qcj)  sin  f3q  cos  mcj)  cos  nd]  d(j), 

^  ^  ./ - 'TT 


q=0 


(101) 
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Q 


C'mn  /  ® 

7 r 


q= o 


cos  (3q  /  cos  cos  mcj)  cos  ncj>d(j> . . . 


■  sin  f3q  /  sin  q<p  cos  mcj)  cos  n<pdcj> 


(102) 

Attention  to  the  integration  limits  reveals  that  the  sin  and  cos  pair  within  the  second 
term  of  (102)  reduce  to  zero,  or 


sin  qcj>  cos  mcj)  cos  rufidtfi  =  0, 


(103) 


which  simplifies  (102)  such  that 


Q 


aq  cos  / 3q  I  cos  qcj)  cos  rruj)  cos  n(j)d(j). 
q= o 


7T 


(104) 


Then,  applying  a  cosine  identity  to  the  m  and  n  indexed  terms  produces 


cos  mcj)  cos  n(j>  =  -  [cos(m  +  n)cj>  +  cos  (m  —  n)<j>\ , 


(105) 


which  can  then  be  substituted  back  into  (104)  as 


Q 


Cmn  ^  ^  ^  ^q  ftq 


q= o 


cos  qcj)  cos  (m  +  n)cj>dcj>  +  /  cos  qcj)  cos  (m  —  n)cj>dcj> 


(106) 


Finally,  the  following  integral  identity  can  be  applied  for  further  reduction, 


cos  mcj)  cos  ncj)dcj)  =  7rSmn , 


(107) 


where  m,  n/0  and  Smn  is  the  Kronecker  delta.  This  leads  to  the  final  computationally 
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efficient  formulation 


where 


1 

2 


Q 

^  ^  Qiq  COS  (3q  [(f^m+n 

q= 0 


“1“  .m— 


q,m—n J  > 


(108) 


{1  . . .  q  —  m  ±  n 

(109) 

0  . . .  q  7^  m  ±  n 


In  expanded  matrix  notation,  this  can  be  practically  implemented  such  that  (for 
m  —  0, . . . ,  M  and  n  —  0, . . . ,  N) 


r< 

w mn 


Coo  C0i  •  •  •  Coat 
Cio  C\  i  . . .  Ciat 


Cmo  Cmi 


C 


MN 


(110) 


This  computationally-efficient  formulation,  relying  on  the  impedance  parameter  ex¬ 
pansion  into  a  Fourier  series,  can  now  be  utilized  during  boundary  condition  enforce¬ 
ment  in  later  sections.  Note,  in  the  case  where  the  impedance  sheet  is  0- invariant, 
the  orthogonal  testing  operator  leads  to  a  simple  single-element  expression,  which 
can  be  computationally  expressed  as  a  vector.  When  a  ^-dependency  is  introduced 
to  the  impedance  sheet,  the  orthogonal  testing  operator’s  ineffectualness  results  in  a 
two-element  expression,  which  can  be  computationally  expressed  as  a  matrix. 


25 


3.2  Free-space  Filled  Tapered  Cylindrical  Impedance  Sheet 

Consider  an  arbitrary  plane  wave  normally  incident  along  the  x  direction  on  a 
cylinder  infinite  in  extent  in  the  5  direction,  as  illustrated  in  Figure  6.  The  cylinder, 
having  radius  a,  is  assumed  to  be  electrically  thin  and  is  comprised  of  ^-dependent 
high-contrast  material  parameters  e, /to ,cr.  Based  on  these  assumptions,  the  cylinder 
can  be  modeled  as  a  generalized  ^-dependent  impedance  sheet.  Given  that  the  dimen¬ 
sions  of  the  physical  scatterer  are  invariant  along  the  z  direction  as  is  the  normally 
incident  plane  wave,  the  total  fields  decouple  into  two  independent  field  sets,  TE2 
and  TMC  The  following  analysis  will  focus  on  the  T1VF  polarization  set,  followed  by 
the  TEZ  case. 


d«  A 


E[  ,H[ 

AAA* 


ES,HS 

*AAA 


E',H' 

AAA* 


G.Mo) 


'  £>  fJ-Q  >  G 


z\  z\ 


+  n  =  p 

zv 

- ►x 

fe»A*o) 
ES,HS 

AAA> 


Figure  6.  Side  and  top  view  of  bistatic  scattering  from  a  cylinder  comprised  of  a 
0-dependent  impedance  sheet. 
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3.2.1  TM~  Polarization  Development. 

The  TMZ  incident  electromagnetic  field  components  expressed  in  computationally 
efficient  form,  as  referenced  from  (77)-(79),  are 

OO 

<t>)  =  E0z  £nj~nJ n(kop)  COS  710,  (111) 

n=0 

Hlp(p,  <t>)  =  E  nr  "Jn(M  sin  710,  (112) 

jkoVoP  “ 

rji  00 

0)  =  —  6nrnJ'n(k0P)  COS  710,  (113) 

{1  if  n  =  0 
2  ifn^O 

where  Po  =  ^\/Po^o,  Vo  =  \Z~PoJt o,  en  is  the  Neumann  number,  Jn  is  the  Bessel 
function  of  the  first  kind  of  order  n  and  H;,2^  is  the  Hankel  function  of  the  second  kind 
of  order  n. 

The  scattered  field  can  be  determined  with  knowledge  of  the  incident  field  and 
the  physical  scatter.  In  this  case,  the  z-invariant  TMZ  set  formulation  will  induce  a 
corresponding  z-directed  surface  current  density  on  the  surface  of  the  impedance  sheet 
cylinder.  This  surface  current  density  will  maintain  the  scattered  field  and  requires 
an  appropriate  mathematical  representation.  With  a  physical  understanding  of  the 
scattering  phenomena,  it  can  be  assumed  that  the  scattered  field  will  naturally  radiate 
outward.  This  necessitates  using  Hankel  functions  of  the  second  kind  to  represent  the 
maintained  outgoing  waves.  Additionally,  utilizing  the  cosine  function  adequately  and 
efficiently  models  the  periodicity  of  the  scattered  field  filling  all  angular  space.  Finally, 
the  scattered  field  is  ultimately  maintained  by  the  incident  field  and  must  therefore  be 
proportional  to  the  incident  field  strength.  Thus,  a  suitable  computationally  efficient 
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representation  to  the  scattered  field  is 


K (p,  0)  =  El0z  E  an£nj  nH^2)  (hop)  cos  n0.  (114) 

n= 0 

The  magnetic  field  components  in  computationally  efficient  form  may  may  then  be 
found  by  substituting  (114)  into  (70),  such  that 

O  T?i  °° 

Hp(p,  0)  =  — —  E  ann3~ntt-n\koP)  sin n(j),  (115) 

jkoVoP  “ 

rri  °° 

Hf{p,  0)  =  —  E  anen3~nft(nY  (kop).  cos  n0.  (116) 

The  physical  cylinder,  being  comprised  of  a  0-' varying  impedance  sheet,  does  not 
scatter  the  entirety  of  the  incident  field.  The  portion  that  is  transmitted  into  the 
cylinder  may  be  adduced  from  the  previously  discussed  z-directed  surface  current 
density  with  the  assumption  that  the  impedance  sheet  is  infinitesimally  thin  (d  <<  A). 
This  transmitted  field  will  become  a  traveling  standing  wave  within  the  cylinder 
and  can  be  appropriately  modeled  with  a  Bessel  function  of  the  first  kind.  The 
same  periodicity  and  proportionality  requirements  still  hold,  resulting  in  the  following 
suitable  computationally  efficient  transmitted  field 


Et(p,  0)  =  E'f)z  E  bn£nj  nJ n(k0p)  COS  77.0.  (117) 

n= 0 

The  magnetic  field  components  in  computationally  efficient  form  may  then  be  found 
by  substituting  (117)  into  (70),  such  that 


2  EL 

jkoVop 


OO 

E  bnnj~nJn(k0p)  sin  n0, 

n=  1 


H‘M  </>) 


(118) 


With  the  electric  and  magnetic  fields  in  all  regions  accounted  for,  the  unknown  ex¬ 
pansion  coefficients  an  and  bn  can  be  found  by  invoking  the  impedance  sheet  boundary 
conditions,  as  referenced  from  (31)  and  (32), 


n  x  E+  —  n  x  E  , 
h  x  (n  x  E±)  =  -Z,(<t>) n  x  (H+  -  H~), 


(120) 

(121) 


where  E±  =  E( a^A),  H±  =  n  is  the  unit  outward  normal,  and  Zs  is  the 

sheet  impedance  which  is  assumed  to  vary  with  cj)  in  this  research.  Substituting  (111), 
(114),  and  (117)  into  (120)  with  h  —  p  results  in 


cnj  nJn(k0a)  cos  ncj)  +  E*0z  ^  anenj  nH^  (k0a)  cos  n<fi 


n= 0 


n= 0 


(122) 


bnenj  nJn(k0a)  cosn0 

ii=0 


Then,  applying  the  testing  operator  (95)  results  in  the  following  (for  m  —  0, . . . ,  oo) 


J2Jn(k0a)^l—  I  cosm(/)cosn(j)d(j)+j2anftin\koa)t-^—  I  cos  m(j)  cos  rufichfi  = 

„  ./ -7 r  „  ./ Tr 


■— n  r  7r 


Cn] 


:—n  /*7r 


n=0 


72=0 


^  ^ ^nJn(^0®) 


;  -n  /*7T 


n= 0 


7T 


cos  m(f)  cos  rupdcj), 


(123) 


where  for  m  =  n  /  0 


cos  rruj)  cos  n(j)d(j)  =  i rSr 


(124) 
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Substituting  (124)  into  (123)  produces 

OO  OO  OO 

E  +  E  =  E  bn^nj  Jn(^O^)^  mri'  (125) 

n= 0  n= 0  n=0 

The  Kronecker  delta  function  effectively  sifts  through  the  infinite  series  summation  to 
identify  only  those  electromagnetic  modes  which  satisfy  both  both  Maxwell’s  equa¬ 
tions  and  the  boundary  conditions,  leading  to  the  following  relation 

Jn(^O^)  +  ^ (^0^)  bn3 n(^k()Cpj .  (126) 

Enforcing  the  second  boundary  condition  will  lead  to  an  additional  relation.  Then, 
with  two  equations  and  two  unknown  expansion  coefficients  an  and  bn,  the  scattered 
field  can  be  determined.  Hence,  substituting  (111),  (113),  (114)  ,  (116),  and  (119) 
into  (121)  with  h  =  p  results  in 

P  x  P  x  4>(El(a,(j))  +Esz(a,<f))) 

-P  x  0  ( El(ai  <t>)  +  0)) 

El(a,(j))  +  Esz(a,(j)) 

which  can  be  further  expanded  as 

°o 

EL  E  W  n  Jn(M  +  anH^(fcoa)  cos n(f>  = 

n=0 

(1  —  bn)J'n(koa )  +  anH^  (koa)  cos  n<f>.  (128) 

The  calculation  of  the  field  expansion  coefficients  an  and  bn  can  be  readily  determined 


=  -ZM)p  x  0  [H;(a,  0)  +  H;(a,  0)  -  H^a,  0)]  , 
=  -Z^Yp  x  0  [H;(a,  0)  +  H}(a,  0)  -  H^a,  0)]  , 
=  Zs(0)  [HA  a,  0)  +  HAa ,  0)  -  HAa,  0)]  ,  (127) 
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by  application  of  the  previously  utilized  testing  operator,  resulting  in 


cos  rruj)  cos  n&dcj)  = 

^  -I  7T  / — 7r 


n= 0 


WO  ^0 


^^enj  n  (1  ~  &„) J^(/coa)  +  anH^  (/c0a)  —  /  Zs(<f))  cos  mcj)  cos  n(j)d(j).  (129) 

^  ^  ./  - TT 


This  can  then  be  further  reduced  by  applying  the  Fourier  series  expansion  of  the 
impedance  sheet,  as  developed  in  (94)-(110),  such  that 


^  ^  ^nj  Jn(^OQ)  T  UnH^(^OQ') 


dinn 


71=0 


Wo  ^0 


y,  enj  n  (1  -  bn)J'n(k0a)  +  anR^y(k0a) 


C„ 


(130) 


With  two  equations,  (126)  and  (130),  and  two  unknown  expansion  coefficients,  an 
and  bn,  the  desired  scattered  field  expansion  coefficient  can  be  readily  found.  Solving 
(126)  for  the  expansion  coefficient  bn  leads  to 


bn  —  IT  cl. 


H,j2)(M 

Jn(^oQ') 


(131) 


which  can  then  be  substituted  back  into  (130)  producing 


Y.i 


J'QoCn] 


J?i(^o^0  anrtn\k0a) 


Vm.n 


71=0 


E 

n= 0 


enrn  |(1-1-  )  J'n(k0a)  T  anR^y(k0a) 

Jn{k0a)  I 


Cmn,  (132) 


y  .T/oTJ  "  Jn(M+«nHi2)(M 


71=0 


&m,n, 


E- 

n= 0 


enJ~ 


Jn(^0®) 


Jn(MH(2)'(M  -  yn(koa)R^(k0a)]  Cmn.  (133) 


This  can  be  additionally  simplified  by  utilizing  Wronskian  and  recurrence  relations, 
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as  developed  in  Appendix  B,  where 


jn(MHi2)'(M  -  j;,(MHi2)(M  =  -^-a  (134) 


This  leads  to  a  simplified  expression  for  the  scattered  field  expansion  coefficient 


OO  OO 

^  ^  Vo^nj  Jn(^0^')^mn  “1“  ^  ^  ^"nVo^nj 

n= 0  n=0 


2Cn 


OO 

E.  _77  ‘Ljy-Jmn 

an^nJ  —j - 777 - y 

7r/coaJn(/c0a) 

n=(J  x  7 


(135) 


Thus,  a  final  formulation  for  the  TM2  scattered  field  expansion  coefficient  can  be 
found  via 


OO 


|  /  v  ttn^nj 
n= 0 


mn 


2  Cr, 


nr)0k0aJn(k0a) 


—  ^  ^  £nj  Jn.(^0®)^?i 


n=0 


(136) 


In  expanded  matrix  notation,  this  can  be  practically  implemented  such  that 


N 


M 


^  ^  Q'nD'mn  ^  ^  Am , 


n=0 


m=0 


Dmn  Cnj 


B%\k0a)5mn  + 


2  C„ 


wr)0k0aJn(k0a) 


Em  —  ^rnj  Jm(^oQ')- 


1 

fa 

o 

o 

fa 

o 

fa 

o 

y 

_ i 

Eq 

D  io  Du  ■  ■  ■  DiN 

CL\ 

= 

Ei 

Dmo  Dmi  ■  ■  ■  Dmn 

Am 

Em 

(137) 

(138) 

(139) 


(140) 


This  result  represents  a  matrix  equation  for  the  unknown  scattered  field  expansion 
coefficients  an  in  a  readily  solvable  Ax  =  b  form.  Note,  truncation  becomes  necessary 
to  some  value  N,  which  is  chosen  in  this  work  to  be  large  enough  to  reach  convergence 
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within  a  tolerance  of  1  x  10~ ' .  With  a  suitable  computationally  efficient  form  achieved, 
the  scattered  field  can  be  determined  for  any  arbitrary  impedance  sheet  distribution 
of  a  cylinder  with  a  given  0  dependence. 

This  formulation  can  be  verified  by  applying  a  special  case.  For  instance,  replacing 
the  ^-dependent  impedance  sheet  with  a  PEC  cylinder  of  identical  dimensions  (Zs  — » 
0),  causes  (136)  to  simplify  to 


^  ^  onenj 

n= 0 


J2^rnUkoa)S 

mm 

n= 0 


(141) 


where  the  scattered  field  expansion  coefficient  an  can  then  be  analytically  solved  for 
as 


Jn(^0®) 

H?M' 


(142) 


This  expression  is  identically  equal  to  the  well-known  TM2  scattered  field  expansion 
coefficient  for  a  PEC  and  grants  more  confidence  for  the  solution. 
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3.2.2  TES  Polarization  Development. 

The  TE"  incident  electromagnetic  field  components  expressed  in  computationally 
efficient  form,  as  referenced  from  (91)-(93),  are 


H'z(p,  0)  =  Hl0z  enj  nJn(k0p)  cos  n0, 

n=0 

(. P ,  0)  =  2^0zVo  V'  nj~nJn(k0p)  Sin  n0, 
jk0p  ^ 


E 


n=  1 


El(p,  0)  =  JHozVo  "Jn( W)  COS  n0. 


— 


n=0 

1  if  n  =  0 

2  if  n  ±  0 


(143) 

(144) 

(145) 


where  ko  =  cu^/po^o,  Vo  =  \J Po/to,  en  is  the  Neumann  number,  Jn  is  the  Bessel 
function  of  the  first  kind  of  order  n  and  H^2'  is  the  Hankel  function  of  the  second  kind 
of  order  n. 

The  scattered  field  can  be  determined  with  knowledge  of  the  incident  field  and 
the  physical  scatter.  In  this  case,  the  z-invariant  TEZ  set  formulation  will  induce 
a  corresponding  z-directed  surface  current  density  on  the  surface  of  the  impedance 
sheet  cylinder.  This  surface  current  density,  though  fictitious,  will  represent  how  the 
scattered  field  is  maintained  and  requires  an  appropriate  mathematical  representation. 
Thus,  a  suitable  computationally  efficient  representation  to  the  scattered  field  is 

OO 

Hz  {p,  0)  =  H0z  anenj~n Hi2)  (hop)  COS  710.  (146) 

n— 0 

The  electric  field  components  in  computationally  efficient  form  may  may  then  be 
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found  by  substituting  (146)  into  (84),  such  that 


Esp(P,<f>) 


25°^°  Y  annj  nH^2)  (kop)  sin  770, 
jko  p 


0)  =  jH^zVo  Y  an6n i  nRnY(kop)  cos 770. 

n=0 


(147) 

(148) 


Similarly,  a  suitable  computationally  efficient  representation  to  the  transmitted  field 
is 


OO 


Ht(p,  0)  =  Hl0z  Y  bntnj  nJn(k0p )  COS  710. 

71=0 

(149) 

The  electric  field  components  in  computationally  efficient  form  may 

found  by  substituting  (149)  into  (84),  such  that 

may  then  be 

£>(P>0)=  • u°z  Ybnni  nJn(W)smri0, 

]k0p  ^ 

(150) 

OO 

El(p,  0)  =  JH0zV 0  X]  6nenj“"Jn(A:op)  cos 710. 

n= 0 

(151) 

With  the  electric  and  magnetic  fields  in  all  regions  accounted  for,  the  unknown  ex¬ 
pansion  coefficients  an  and  bn  can  be  found  by  invoking  the  impedance  sheet  boundary 

conditions,  as  referenced  from  (31)  and  (32), 

n  x  E+  =  n  x  E  , 

fix  (fix  £*)  =  —Zs((f))n  x  (H+  -  H~), 

(152) 

(153) 

where  =  E(a±,cj)),  H ±  =  H(a±,cj)),  n  is  the  unit  outward  normal,  and  Zs  is  the 
sheet  impedance  which  is  assumed  to  vary  with  0  in  this  work.  Substituting  (144), 
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(145),  (147),  (148),  (150),  and  (151)  into  (152)  with  h  —  p  results  in 


nj  nJ'n(^0 p)  COS  710  +  j H^Jj  nH^  (/c0p)  cos  ncj) 


n= 0 


n=0 


mLvY,  b„e„j  nl'n(kap)  cosn<t> 


■  (154) 


71—0 


The,  applying  the  testing  operator  (95)  results  in  the  following  (for  m  —  0, . . . ,  oo) 


5ZJ n(M^y—  /  cosm0cosn0d0  +  ^anH(2)'(/cop)^^  /  cos  m<j)  cos  n<j)d(f)  = 

^  — 7T  ^  — 7T 


;  —  n  rn 


;  —  n  /*7r 


71—0 


where  for  m  =  n  ^  0 


n=0 


&™Jn(W) 


CnJ 


n  /■7T 


n=0 


7T 


cos  mcj)  cos  rufidtfi 


(155) 


cos  mcj)  cos  n(f)d(f)  =  Trip 


(156) 


Substituting  (156)  into  (155)  produces  the  following  relation 


yn(k<>a)  +  a„H< fy(k0a)  =  6ny(A0o). 


(157) 


Applying  the  second  boundary  condition  will  lead  to  an  additional  relation.  Then, 
with  two  equations  and  two  unknown  expansion  coefficients  an  and  bn,  the  scattered 
field  can  be  determined.  Hence,  substituting  (143)-(145),  (146)-(148),  and  (149)-(151) 
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into  (153)  with  n  =  p  results  in 


P  x 


p  x  ( pEp(a ,  0 )  +  (pE^a,  0) ) 
p  x  i  (EJ(a,0)  +  ^(a,0)) 


M(a,0)  +  £i(a,0) 


-Zs(0)p  x  5  [if*(a,  0)  +  iff (a,  0)  -  iff (a,  0)]  , 

(158) 

-Zs(0)p  x  5  [Hl(a,  0)  +  iff (a,  0)  -  iff  (a,  0)]  , 

(159) 

-Zs(0)  [iff (a,  0)  +  Hsz(a,  0)  -  if* (a,  0)]  , 

(160) 


which  can  be  further  expanded  as 


JH0zV^2enj  n  J'n(k0a)  +  anR(*y (k0a)  cos (n0)  = 

00 

H^zZs(<P)J2^rn  [(1  -  &n)J»(M  +  anH^(fc0a) 


n=0 


n=0 


cos  (n0).  (161) 


Applying  the  testing  operator  (95)  and  the  Fourier  series  expansion  of  the  impedance 
sheet  (96),  results  in 


JV^nJ 


n= 0 


Jn(^o^)  ^  (koa) 


&m.n, 


^  ^  (1  ^n)Jn(^O^)  “1“  ^0^) 


n=0 


n 

wmn* 


(162) 


With  two  equations,  (157)  and  (162),  and  two  unknown  expansion  coefficients,  an 
and  bn,  the  desired  scattered  field  expansion  coefficient  can  be  readily  found.  Solving 
(157)  for  the  expansion  coefficient  bn  leads  to 


b 


n 


1  +  CLn 


Hl2)/(M 

Jn(^0®) 


(163) 
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which  can  then  be  substituted  back  into  (162)  producing 


Y.J 


JVo^nJ 


J  n(koa)  +  anR^y(k0a) 


3 m.n . 


71—0 


E 

n= 0 


enJ 


1  —  1  —  a 


H  W(k0a) 

Jn(^OQ') 


71  -r/ 


Jn(/c0a)  +  anH^2)(fc0a) 


a 


(164) 


y,  jyoCnj  n  yn(k0a)  +  anR^y  (k0a) 


71—0 


&m,n, 


y  ttn  y 


enJ~ 


71—0 


^(M  L 


-  Jn(MHi2)'(M 


Cn 


(165) 


This  can  be  additionally  simplified  by  utilizing  Wronskian  and  recurrence  relations 
developed  in  Appendix  B,  where 


yn(k0a)R^\k0a)  -  Jn(MHi2)'(M  =  •  (166) 


This  leads  to  a  simplified  expression  for  the  scattered  field  expansion  coefficient 


OO  OO 

-  y  voenj-nyn  (k0°)smn  -  y  ^nVo^nj  Hn  ^  (^0®)^  77171 
71—0  71—0 


E 

71—0 


^n^nj 


2Cn 


7r/coaJ^(/coa) 


(167) 


Thus,  a  final  formulation  for  the  TE~  scattered  field  expansion  coefficient  can  be 
found  via 


OO 


/  v  ttn^n] 


71—0 


H  W(k0a)5mn  + 


2  C„ 


nri0k0ayn(koa)  \ 


OO 

y  enj~r\l,n(k()a)Sv 

n= 0 


(168) 


In  expanded  matrix  notation,  this  can  be  practically  implemented  in  (137)  such  that 
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Dmn  Cnj 


( koa)5mn  + 


2  C„ 


irr}okoaJ'n(k0a) 


Em,  £ mj  ■ 


(169) 

(170) 


This  formulation  can  be  verified  by  applying  a  special  case.  For  instance,  replacing 
the  ^-dependent  impedance  sheet  with  a  PEC  cylinder  of  identical  dimensions  (Zs  — > 
0),  causes  (168)  to  simplify  to 


E 

71—0 


H  W(k0a)6„ 


-E 


W" 


lJ'n(k0a)Sn 


(171) 


n= 0 


where  the  scattered  field  expansion  coefficient  an  can  then  be  analytically  solved  for 
as 


Qjn 


J'(M 


(172) 


This  expression  is  identically  equal  to  the  well-known  TEZ  scattered  held  expansion 
coefficient  for  a  PEC  and  grants  more  confidence  for  the  solution. 
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3.3  PEC  Filled  Tapered  Cylindrical  Impedance  Sheet 


Consider  a  plane  wave  traveling  along  the  forward  x  direction  normally  incident  on 
a  cylindrical  structure  that  is  infinite  in  extent  along  z,  as  illustrated  in  Fig.  7.  The 
structure  is  comprised  of  an  electrically-thin  (d  <C  A)  cylindrical  impedance  sheet 
having  radius  cq  and  an  inner  solid  pedestal  of  radius  a2.  The  ^-dependent  high- 
contrast  material  parameters  of  the  impedance  sheet  are  e,po,cr  and  the  pedestal  is 
assumed  to  be  comprised  of  a  PEC.  The  regions  a2  <  p  <  «i  and  p  >  cq  are  assumed 
to  be  free-space  filled. 


/v 

7 


Figure  7.  Side  and  top  view  of  bistatic  scattering  from  a  cylinder  comprised  of  a 
0-dependent  impedance  sheet  wrapped  about  a  PEC  pedestal. 

Since  both  the  scattering  structure  and  incident  excitation  are  ^-invariant,  Maxwell’s 
equations  decouple  into  two  independent  field  sets,  TEZ  and  TMC  In  addition,  based 
on  the  above  assumptions,  an  impedance  sheet  boundary  condition  may  be  utilized 
to  accommodate  the  scattering  effects  of  the  impedance  sheet.  The  z-invariance  and 
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impedance  sheet  boundary  condition  implementation  greatly  facilitates  the  analytical 
development  for  the  bistatic  scattering  width  calculations,  as  discussed  next. 


3.3.1  TM2  Polarization  Development. 

For  TM2  polarization,  the  tangential  electric  and  magnetic  held  (i.e. ,  incident  plus 
scattered)  components  in  the  region  p  >  a\  can  be  expressed  according  to  (77)-(79) 
and  (114)-(116)  as  (assuming  an  time  dependence) 


ET (p,  0)  =  E'0z  enJ~n  [-UM  +  anliW(k0p) 


Ek 


71—0 

oo 


K'(p,  </>)  =  E  ' '(m 


™  ^ 0 


cos  n0, 
cos  n0, 


(173) 

(174) 


where  k0  =  LUy/e0p0,  r]0  =  \J po/eo,  en  is  the  Neumann  number,  Jn  is  the  Bessel 
function  of  the  first  kind  of  order  n,  is  the  Hankel  function  of  the  second  kind 
of  order  n  and  the  prime  represents  differentiation  with  respect  to  the  argument. 
Similarly,  the  tangential  fields  in  regions  <22  <  p  <  a\  may  be  written  as 


Et(p,  0)  =  E'0z  ^  enj  n  bnJn(k0p )  +  cnR^\k0p) 


Ek 


71—0 

OO 


<t>)  =  —  y,  <*T"  Wi(M  +  CnH'P’M 

3  Vo  ^ 


71—0 


cos  n0, 
cos  n0. 


(175) 

(176) 


The  unknown  amplitude  coefficients  an ,  bn ,  and  cn  are  found  via  boundary  condi¬ 
tion  enforcement  at  p  =  a\  and  p  =  a2.  Note,  the  expansion  coefficient  of  particular 
relevance  for  bistatic  scattering  width  calculations  is  an.  Enforcing  continuity  of  the 
tangential  electric  and  magnetic  fields  at  the  PEC  interface  (i.e.,  p  =  a-2)  and  applying 
the  suitable  testing  operator  results  in 


bn^n{koCi2)  +  CnH^(/c0a2)  —  0. 


(177) 
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Next,  it  is  necessary  to  enforce  the  impedance  sheet  boundary  conditions 


n  x  E+  =  n  x  E  , 


(178) 


ft  x  (ft  x  J±)  =  -Z,(<f)h  x  (H+  -  H~), 


(179) 


where  =  E(af,<l>),  H±  =  h  is  the  unit  outward  normal,  and  Zs  is  the 

sheet  impedance  which  is  assumed  to  vary  with  0  in  this  work.  Substituting  (173)- 
(174)  and  (175)-(176)  into  (178)  and  (179)  with  n  =  p  and  applying  the  suitable 
testing  operator  leads  to 


Jn(Ml)  +  an^-n\^  0al)  —  ^nJn(^Oal)  +  C„H^  (&0al) , 


(180) 


y, :imenj  n  Jn(^0ai)  +  anH^(fc0ai) 


71—0 


&m.n. 


y  enj  n  (1  -  &n)JU^0ai)  +  anH^^fcoai)  -  cnR%y (ko^) 


n= 0 


a 


(181) 


where  the  substitution  of  the  Fourier  series  expansion  of  the  0-varying  impedance 
sheet  has  been  included.  Solving  (177)  in  terms  of  cn  results  in 


Cn  hr 


Jn(^0°2) 

H  l?(k0a2y 


(182) 


which  can  then  be  substituted  into  (180)  leading  to  the  expression 


bn  = 


Jn(Ml)  +  (&0al) 


(183) 


The  desired  TMZ  scattered  held  expansion  coefficients  an  are  readily  solved  by 
substituting  (182)  and  (183)  into  (181),  leading  with  the  inclusion  of  Wronskian  and 
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recurrence  relations  to  (for  m  —  0, . . . ,  N ) 


n= 0 
oo 

^  ^  cinenj 


7T  ^O'^O a  1  Jn  (  ^0 a  1 )  <5mn  T 


2Jn(^0®2)C*r7 


Jn(^oai)H^2)(/c0a2)  -  Jn(k0a2)Hi;Z\k0al) 


n— 0 


7t^o^0ohH^2)  (k0ai)Smn  + 


2H  W(k0a2)C„ 


Jn(k0a1)H(n2)(k0a2)  -  Jn(k0a2)H(n2)(k0a1) 


(184) 

In  expanded  matrix  notation,  this  can  be  practically  implemented  as  in  (137).  This 
formulation  can  be  verified  analytically  by  employing  a  special  case.  If  the  PEC 
tube  within  the  cylinder  comprised  of  a  ^-dependent  impedance  sheet  is  removed, 
then  (184)  should  identically  reduce  to  (136).  This  can  be  readily  accomplished  by 
allowing  the  imposed  PEC  boundary  condition,  as  referenced  by  (182),  to  equal  zero, 
which  leads  to  the  following 


OO 

n= 0 


Hi2\k0a1)dmn  + 


2  a 


nr}ok0aiJn(koai) 


OO 

^  ^  ^nj  *1  n ( /cqu  ] ) Srnn .  (185) 

n= 0 


Thus,  as  the  PEC  tube  within  the  cylinder  comprised  of  a  ^-dependent  impedance 
sheet  is  removed,  then  (184)  reduces  identically  to  (136)  and  grants  more  confidence 
for  the  solution. 


3.3.2  TES  Polarization  Development. 

For  TEZ  polarization,  the  tangential  electric  and  magnetic  held  (i.e. ,  incident  plus 
scattered)  components  in  the  region  p  >  cq  can  be  expressed  according  to  (91)-(93) 
and  (146)-(148)  as  (assuming  an  time  dependence) 


h^(pa)  = 

n= 0 
oo 

E‘p(P,<b)  =  3K,nY.^" 

71=0 


J n(k0p)  +  anenj  "H )2)(/c0p) 

J„(M  T  cin6nj  ^  (k0p) 


cos  n0, 
cos  ncj). 


(186) 

(187) 
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where  k0  =  LU^/e0po,  rj0  =  \Zpo/e0,  en  is  the  Neumann  number,  Jn  is  the  Bessel 
function  of  the  first  kind  of  order  n,  is  the  Hankel  function  of  the  second  kind 
of  order  n  and  the  prime  represents  differentiation  with  respect  to  the  argument. 
Similarly,  the  tangential  fields  in  regions  <22  <  P  <  a\  may  be  written  as 

00 

Hl{p,  4>)  —  Hqz  enj  n  bn3n(k0p) +  cnR^\k0p)  cos  n0,  (188) 

n= 0 

00 

E^p,  0)  =  jHl0zp 0  enj~n  bnyn(k0p)  +  cnH^2)'  (k0p)  cos  n(j).  (189) 

)i=0 

The  unknown  amplitude  coefficients  an,  bn,  and  cn  are  found  via  boundary  condi¬ 
tion  enforcement  at  p  =  ci\  and  p  =  0,2-  Note,  the  expansion  coefficient  of  particular 
relevance  for  bistatic  scattering  width  calculations  is  an.  Enforcing  continuity  of  the 
tangential  electric  and  magnetic  fields  at  the  PEC  interface  (i.e.,  p  =  <12)  and  applying 
the  suitable  testing  operator  results  in 

bnJn(koa2)  +  cnH ^  (k0a2)  =  0.  (190) 

Next,  it  is  necessary  to  enforce  the  impedance  sheet  boundary  conditions 

n  x  E+  —  h  x  E~ ,  (191) 

n  x  (n  x  E±)  =  -Z8{<j>)h  x  ( H+  -  H~),  (192) 

where  =  E(af,(f>),  H ±  =  H(af  ,(/>).,  h  is  the  unit  outward  normal,  and  Zs  is  the 
sheet  impedance  which  is  assumed  to  vary  with  <fi  in  this  work.  Substituting  (186)- 
(187)  and  (188)-(189)  into  (191)  and  (192)  with  h  =  p  and  applying  the  suitable 


44 


testing  operator  leads  to 


Jn(^Oal)  +  flnH®  (Ml)  —  ^reJn(^0°l)  +  (&0al)> 


(193) 


Ei 


JVo^nJ 


Jn(^oai)  +  onH(l2')  (A;0ai) 


$mn 


n= 0 


E  e«7  n  (1  -  MJn(Mi)  +  anU.^\koai)  -  cnH^2)(Mi) 


n=0 


a, 


(194) 


where  the  substitution  of  the  Fourier  series  expansion  of  the  0-varying  impedance 
sheet  has  been  included.  Solving  (190)  in  terms  of  cn  results  in 


J  'n(k0a2) 

"Hf'lW 


(195) 


which  can  then  be  substituted  into  (193)  leading  to  the  expression 


h;  — 


Jn(Ml)  ~ 


(196) 


The  desired  TMZ  scattered  held  expansion  coefficients  an  are  readily  solved  by 
substituting  (195)  and  (196)  into  (194),  leading  with  the  inclusion  of  Wronskian  and 
recurrence  relations  to  (for  m  —  0, . . . ,  N ) 


E 

n= 0 


71  hyrioci \ ^n(k()(ii ^Smn  T 


2Jn(^0a2)Cr) 


^  ^ an^nj 


n= 0 


7rMoaiH^2)'(Mi)(5mn  +  — 


2Hi2),(M2)C'mn 


(197) 

In  expanded  matrix  notation,  this  can  be  practically  implemented  as  in  (137).  This 
formulation  can  be  verified  analytically  by  employing  a  special  case.  If  the  PEC 


45 


tube  within  the  cylinder  comprised  of  a  ^-dependent  impedance  sheet  is  removed, 
then  (197)  should  identically  reduce  to  (168).  This  can  be  readily  accomplished  by 
allowing  the  imposed  PEC  boundary  condition,  as  referenced  by  (195),  to  equal  zero, 
which  leads  to  the  following 


E 

77—0 


(fcoal)<5mn  + 


2  C„ 


7T77o^OalJn(^Oal) 


y,  enj  nyn(koai)5mn.  (198) 

77—0 


Thus,  as  the  PEC  tube  within  the  cylinder  comprised  of  a  ^-dependent  impedance 
sheet  is  removed,  then  (197)  reduces  identically  to  (168)  and  grants  more  confidence 
for  the  solution. 
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3.4  Dielectric  Filled  Tapered  Cylindrical  Impedance  Sheet 

Consider  a  plane  wave  traveling  along  the  forward  x  direction  normally  incident 
on  a  cylindrical  structure  that  is  infinite  in  extent  along  z,  as  illustrated  in  Fig. 
8.  The  structure  is  comprised  of  an  electrically-thin  (d  <C  A)  cylindrical  impedance 
sheet  having  radius  cq  and  an  inner  solid  pedestal  of  radius  a2.  The  ^-dependent 
high-contrast  material  parameters  of  the  impedance  sheet  are  e,  /i,  a  and  the  pedestal 
is  assumed  to  be  comprised  of  simple  (i.e.,  linear,  homogeneous,  isotropic)  medium 
€d,  Hd-  The  regions  a2  <  p  <  a\  and  p  >  a±  are  assumed  to  be  free-space  filled. 

/v 
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Figure  8.  Side  and  top  view  of  bistatic  scattering  from  a  cylinder  comprised  of  a 
^-dependent  impedance  sheet  wrapped  about  a  dielectric  pedestal. 

Since  both  the  scattering  structure  and  incident  excitation  are  ^-invariant,  Maxwell’s 
equations  decouple  into  two  independent  field  sets,  TEZ  and  TMT  In  addition,  based 
on  the  above  assumptions,  an  impedance  sheet  boundary  condition  may  be  utilized 
to  accommodate  the  scattering  effects  of  the  impedance  sheet.  The  2- invariance  and 
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impedance  sheet  boundary  condition  implementation  greatly  facilitates  the  analytical 
development  for  the  bistatic  scattering  width  calculations,  as  discussed  next. 


3.4.1  TM2  Polarization  Development. 

For  TM2  polarization,  the  tangential  electric  and  magnetic  held  (i.e. ,  incident  plus 
scattered)  components  in  the  region  p  >  a\  can  be  expressed  as  (assuming  an  e^ut 
time  dependence) 


Et’s(p,  0)  =  Eiz  Y  enj-n  [j„(M  +  anHf  (M 


Ek 


n= 0 
oo 


K'(p,  ‘P)  =  —  Y.  ‘"j- "  J»(M  +  m 


™  ^ 0 


cos  nd, 
cos  nq \ 


(199) 

(200) 


where  k0  =  LUy/e0p0,  p0  =  \J /i0 /e0 ,  en  is  the  Neumann  number,  Jn  is  the  Bessel 
function  of  the  first  kind  of  order  n,  is  the  Hankel  function  of  the  second  kind 
of  order  n  and  the  prime  represents  differentiation  with  respect  to  the  argument. 
Similarly,  the  tangential  fields  in  regions  a-2  <  p  <  a\  and  p  <  a,2,  respectively,  may 
be  written  as 


El(p,  </>)  =  EL  X  E»7~" 

n=0 


hI(p,4> ) 


OO 


n= 0 


bnJn(kop)  +  cnHn\k0p) 

KJ'nikop)  +  cnR^y(k0p) 


cos  nq i, 

1  cos  n(j), 


OO 


(201) 

(202) 


Ez(p ,  0)  =  El0z  Y  dntnj  nJ n(kdp)  cos n0,;  (203) 

n=0 
T?i  °° 

H}(p,  0)  =  —  Y  dn£nj~nJ'n(kdp)  COS  710,  (204) 

where  kd  =  u^edpd  and  rjd  =  \J pdjed. 

The  unknown  amplitude  coefficients  an,  bn,  cn,  and  dn  are  found  via  boundary 
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condition  enforcement  at  p  —  a\  and  p  =  a 2.  Note,  the  expansion  coefficient  of  par¬ 
ticular  relevance  for  bistatic  scattering  width  calculations  is  an.  Enforcing  continuity 
of  the  tangential  electric  and  magnetic  fields  at  the  pedestal  interface  (i.e.,  p  =  a2) 
results  in 


EL 


E 

n=0 


tnf 


&«Jn(^0a2)  +  CrMn\^0a2) 


COS  ncf) 


OO 

ELY,  dnenj  nJn{kda2)  cos n<f>, 

n= 0 


(205) 


Eqz 

JVo 


J2br 

n= 0 


tnf 


/  1  E i  A 

bn3'n(k0a2)  +  cnR^y(k0a2)  cos ruf)  =  — —  dnenj~nJ'n(kda2 )  cos ncf), 

J  dVd  t'o 


(206) 


and  applying  the  suitable  testing  operator  to  (205)  and  (206)  leads  to  the  following 
(for  n  =  0, ,  cx)) 


Vd 


’bn^n(koCL2)  "h  CnH()(fc0a2)  dn Jn {kdCL‘2) , 

bnyn(k0a2)  +  cnR^y  (k0a2)]  =  Podn3'n{kda2) . 


(207) 

(208) 


Solving  (207)  and  (208)  in  terms  of  dn  and  equating  leads  to  the  expression 


Cn  bn 


P0^n(k()Cl2) ^ n{kd&2)  Pd^ n{koCl2)^n{kdCJ/2) 
:p0Hn\k0a2)yn{kda2)  -  pdR{p;y  (koa2)Jn(kda2) 


(209) 


Next,  it  is  necessary  to  enforce  the  impedance  sheet  boundary  conditions 


fi  x  E+  =  n  x  E  , 
fix  (fix  E*)  =  —Zs(4>)n  x  (. H+  -  H~), 


(210) 

(2H) 


where  E ±  =  _E(af,0),  H*1  =  H(af,<f>),  h  is  the  unit  outward  normal,  and  Zs  is  the 
sheet  impedance  which  is  assumed  to  vary  with  cj)  in  this  work.  Substituting  (199)- 
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(200)  and  (201)-(202)  into  (210)  and  (211)  with  h  =  p  and  applying  the  suitable 
testing  operator  leads  to 


Jn(M)Oi)  +  anR^\k0ai)  =  bnJn{kvai)  +  cnB^\k0ai),  (212) 

OO 

^  ^  jVO^nj  Jn(^O^l)  “1“  ^ O^l)  $mn 

«,  •  (213) 

jAnV”  (1  -  f>„)y(*0«l)  +Or.H®'(*0°l)  -CnHy'fMl)  Cm„ 

n= 0 

An  expression  for  the  expansion  coefficient  bn  in  terms  of  an  may  be  found  by  substi¬ 
tuting  (209)  into  (212)  and  solving  to  produce 


Jn(^0®l) 


Jn  (^0®l)  “h  «nHi  ^(&0®l) 

hO  Jn(^0®2)  Jn(^d®2)  Vd^ni^O^'2)  Jrc(^ci®2)  jj-(2)  /  r,  \ 


•  (214) 


The  desired  TMZ  scattered  held  expansion  coefficients  an  are  readily  solved  by 
substituting  (209)  and  (214)  into  (213),  leading  to  (for  m  —  0, . . . ,  N) 


(215) 

In  expanded  matrix  notation,  this  can  be  practically  implemented  as  in  (137).  Also 
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note  that  the  special  cases  of  a  cylindrical  impedance  sheet  filled  with  free-spaee  or  a 
cylindrical  PEC  are  easily  obtained  by  letting  r]d  =  oo  or  77^  =  0,  respectively. 


3.4.2  TE2  Polarization  Development. 

For  TM2  polarization,  the  tangential  electric  and  magnetic  held  (i.e.,  incident  plus 
scattered)  components  in  the  region  p  >  a\  can  be  expressed  as  (assuming  an  e^ut 
time  dependence) 


0)  =  H'0z  tnJ~n  [j„(M  +  anR^\k0p) 


COS  710, 


n= 0 


E^(p,  4)  =  jH0zV0  e 'nj  n  Jn(^0 P)  +  OnH ^  (k0p) 


COS  710, 


71=0 


(216) 

(217) 


where  k0  =  LUy/e0p0,  p0  =  \J /i0 /e0 ,  en  is  the  Neumann  number,  Jn  is  the  Bessel 
function  of  the  first  kind  of  order  n,  is  the  Hankel  function  of  the  second  kind 
of  order  n  and  the  prime  represents  differentiation  with  respect  to  the  argument. 
Similarly,  the  tangential  fields  in  regions  a-2  <  p  <  a\  and  p  <  0,2,  respectively,  may 
be  written  as 


Hl(p,  0)  =  H^^enj  n  bn3n(k0p)  +  cnR^\k0p)  cos  n0, 

71=0 

OO 

Kip,  0)  =  jHbzV 0  Y  en^  n  b>Mkop)  +  cjl^'ikop)  COS  710. 

71=0 

OO 

Hdz(p ,  0)  =  Hqz  Y  €njndnJn(kdp)  COS  77.0, 

71=0 

OO 

E}(p,  0)  =  jHq. zVd  Y  enj~ndnyn{kdp)  cos  710. 


n= 0 


(218) 

(219) 

(220) 
(221) 


where  kd  =  u^edpd  and  rjd  =  s]pd/ed. 

The  unknown  amplitude  coefficients  an,  bn,  cn,  and  dn  are  found  via  boundary 
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condition  enforcement  at  p  —  a\  and  p  =  a 2.  Note,  the  expansion  coefficient  of  par¬ 
ticular  relevance  for  bistatic  scattering  width  calculations  is  an.  Enforcing  continuity 
of  the  tangential  electric  and  magnetic  fields  at  the  pedestal  interface  (i.e.,  p  =  a2) 
results  in 


■iHozrh)  X/  erd  n  bnJ'n(k0a2)  +  cnH,(t2)'  (k0a2) 


n= 0 


cos  ruj)  = 


JHoz>ld  E  enj  ndnJ'n(kda2)  cos  mp 


(222) 


n= 0 


H()z  ^  ^  ^ nj  ^nJn(^0^2)  “1“  ^(^0^2) 


n=0 


cos  rub  = 


EE£  nj  ndnJn(kda2)  cosnj) 


(223) 


n= 0 


and  applying  the  suitable  testing  operator  to  (222)  and  (223)  leads  to  the  following 
(for  11  00) 


Vo 


bnyn{k0a2 )  +  cnR^y(k0a2)  =  VddnJ'n(kda2) , 
bn^n^kott^  "h  CnR^ik^a^,)  dn^n{kdci2). 


(224) 

(225) 


Solving  (224)  and  (225)  in  terms  of  dn  and  equating  leads  to  the  expression 


C'n  bn 


Vo^  ni.ko^')  Jn(,kda2)  VdUko02)Ukd02) 

I]0n^y (k0a2)Jn(kda2)  -  r)dH.£)(k0a2)J,n{kda2) 


Next,  it  is  necessary  to  enforce  the  impedance  sheet  boundary  conditions 


(226) 


h  x  E+  =  n  x  E  , 
h  x  (n  x  E±)  =  -Zs(<j>)n  x  ( H+  -  H~), 


(227) 

(228) 


where  EE  =  E(af,(j)),  H±  =  H(af,q !>),  n  is  the  unit  outward  normal,  and  Zs  is  the 


52 


sheet  impedance  which  is  assumed  to  vary  with  0  in  this  work.  Substituting  (216)- 
(217)  and  (218)-(219)  into  (227)  and  (228)  with  h  —  p  and  applying  the  suitable 
testing  operator  leads  to 


yn(k0ai)  +  anB^'ik^)  =  bnJ'n(k0ai)  +  cnR^y  (k0ai) ,  (229) 

OO 

^  ^  jVo^nj  Jn(fcbai)  +  cin ^  (fc0ai)  Smn 

co  •  (230) 

^  ^  ^ nj  (1  ^n)Jri(^0®l)  d"  Cln (/uqO^)  Cn (/ugfl^)  Cmn 

n= 0 

An  expression  for  the  expansion  coefficient  bn  in  terms  of  an  may  be  found  by  substi¬ 
tuting  (226)  into  (229)  and  solving  to  produce 


Jn(*bai ) 


Jn(^oai)  +  onH (k0ai) 

^  n{kdOJ2)  n{k^2)  J«(^rf^2)  ^  \ 

77oH^,(^o«2)Jn(^d«2)  —  VdJin\koa2)J'n(kda2)\  ”  01 


(231) 


The  desired  TMZ  scattered  held  expansion  coefficients  an  are  readily  solved  by 
substituting  (226)  and  (231)  into  (230),  leading  to  (for  m  —  0, . . . ,  N) 


(232) 
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In  expanded  matrix  notation,  this  can  be  practically  implemented  as  in  (137).  Also 
note  that  the  special  cases  of  a  cylindrical  impedance  sheet  filled  with  free-spaee  or  a 
cylindrical  PEC  are  easily  obtained  by  letting  r]d  =  oo  or  rjd  =  0,  respectively. 


54 


IV.  Experimental  Results 


In  the  previous  chapter,  the  scattered  held  expansion  coefficient  development  of 
TM2  and  TE2  polarizations  for  a  cylinder  comprised  of  a  ^-dependent  impedance 
sheet  were  derived.  This  relied  on  general  impedance  sheet  boundary  conditions  in 
addition  to  a  Fourier  series  representation  of  the  impedance  sheet  distribution,  with 
the  utilization  of  Wronskian  and  recurrence  relations,  to  produce  a  computationally 
efficient  matrix  solution  for  the  scattered  held  expansion  coefficients  an.  The  scat¬ 
tering  coefficients  an  can  now  be  included  in  the  expression  for  a  two-dimensional 
scattering  width,  allowing  for  the  scattering  patterns  of  arbitrary  diameter  cylinder 
impedance  sheet  distributions.  The  ability  to  adapt  sheet  impedances  to  effectively 
control  the  overall  nulls  in  the  bistatic  scattering  pattern  will  be  demonstrated,  with 
practically  implemented  impedance  characteristics.  Lastly,  a  novel  rectangular  X- 
band  waveguide  system  will  be  utilized  to  characterize  the  electromagnetic  properties 
of  the  ^-dependent  impedance  sheet. 


4.1  Bistatic  Scattering 


In  order  to  obtain  a  meaningful  representation  of  the  scattered  held,  the  two- 
dimensional  scattering  width  will  be  utilized.  This  will  allow  for  insight  on  how  the 
impedance  distribution  about  the  cylinder  affects  the  scattering  pattern.  The  bistatic 
scattering  width  may  be  computed  via 


C-2D 


lim 

p — >oo 


£s(p,0)r 
£*(p,0)|2J  ’ 


(233) 
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which  upon  using  the  asymptotic  expansion  for  the  scattered  held  becomes  (normal¬ 
ized  with  respect  to  wavelength) 

2 

(234) 

As  previously  mentioned,  the  number  of  terms  N  in  (234)  is  chosen  such  that  con¬ 
vergence  is  reached  to  within  a  tolerance  of  1  x  10“'.  Calculation  of  the  expansion 
coefficients  an  and  corresponding  bistatic  scattering  width  for  various  impedance  sheet 
distributions  is  demonstrated  next. 

In  order  to  demonstrate  the  effect  that  impedance  sheet  tapering  has  on  the 
bistatic  scattering  pattern  (when  compared  to  a  uniform  sheet),  two  scattering  sce¬ 
narios  are  investigated  here.  The  two  cases  involve  a  tapered  impedance  sheet  filled 
with  a  dielectric  and  PEC  pedestal,  respectively.  Note,  as  discussed  in  the  theory,  an 
rc-clirected  normally-incident  plane  wave  excitation  is  assumed. 


&2D 


2 
7 r 


N 


enan  cos  n(f> 


71= 0 


4.1.1  Dielectric  pedestal  results. 

Consider  a  non- magnetic  polystyrene  pedestal  having  permittivity  =  1.5eo  and 
radius  a±  =  A  surface- wrapped  (i.e.,  a,2  =  a\)  with  the  tapered  impedance  sheet  dis¬ 
tribution  Zs{(p)  shown  in  Figure  9.  This  distribution  represents  a  physically-realizable 
sinusoidal  impedance  taper  from  0  to  1000-j200  consisting  of  a  single  period  about 
the  cylinder.  As  will  be  shown  below,  this  relatively  simple  complex  impedance  sheet 
distribution  can  effectively  introduce  specihc  nulls  in  the  bistatic  scattering  width 
unattainable  with  a  constant  impedance  sheet  and  allows  the  null  positions  to  be 
readily  controlled. 

The  bistatic  scattering  widths  illustrated  in  Figure  10  depict  the  effects  of  the 
TM*  polarization.  At  its  initial  setup  state,  the  single  sinusoid  impedance  sheet 
distribution  introduces  particular  nulls  near  60  and  100  degrees.  These  nulls  are  non- 
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Figure  9.  Single  sinusoid  (^-dependent  impedance  sheet  distribution. 

trivial,  with  a  two-dimensional  scattering  width  improvement  of  about  17  dB  and 
14  dB,  respectively.  Rotating  the  pedestal  by  23.4  degrees  deepens  the  null  at  60 
degrees  to  a  reduction  of  37  dB  while  maintaining  the  100  degree  null  at  about  12 
dB.  Then,  rotating  the  pedestal  by  90  degrees  effectively  repositions  the  two  nulls, 
resulting  in  two-dimensional  scattering  width  improvements  at  25,  70,  110,  and  150 
degrees  with  respective  reductions  of  about  10  dB,  9.5  dB,  9  dB,  and  7.5  dB.  Overall, 
the  bistatic  scattering  widths  presented  Figure  10  demonstrate  three  key  features  of 
tapered  over  uniform  impedance  sheet.  It  shows  how  a  tapered  impedance  sheet  not 
only  introduces  controllable  and  particular  nulls  in  the  bistatic  scattering  widths,  but 
the  ability  to  significantly  increase  the  reduction  amount  at  certain  angles.  These 
effects  will  be  seen  again  in  the  remaining  analysis. 

Figure  11  illustrates  the  bistatic  scattering  widths  of  TE2  polarization.  At  its 
initial  setup  state,  the  single  sinusoid  impedance  sheet  distribution  introduces  a  broad 
null  centered  at  60  degrees  and  a  particular  null  near  120  degrees.  These  nulls  are 
non-incidental,  with  the  broad  null  reducing  the  two-dimensional  scattering  width  by 
at  least  10  dB  from  about  45  to  75  degrees,  with  a  maximum  reduction  of  18  dB. 
Meanwhile,  the  particular  null  at  120  degrees  improves  the  bistatic  scattering  width 
by  16  dB.  Then,  rotating  the  pedestal  by  90  degrees  effectively  repositions  the  nulls 
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Figure  10.  2D  Bistatic  Scattering  Width  of  TM2  polarization  (Top:  0  degrees,  Middle: 
23.4  degrees,  Bottom:  90  degrees)  for  the  (/(-dependent  impedance  sheet  filled  with  a 
polystyrene  pedestal. 
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Figure  11.  2D  Bistatic  Scattering  Width  of  TEZ  polarization  (Top:  0  degrees,  Middle: 
90  degrees,  Bottom:  104.4  degrees)  for  the  ^dependent  impedance  sheet  filled  with  a 
polystyrene  pedestal. 
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to  about  90  and  135  degrees,  with  corresponding  reductions  of  24  dB  and  21  dB. 
Additionally,  at  this  rotation  a  14  dB  null  remains  at  50  degrees  while  a  9  dB  null 
at  170  degrees  is  introduced.  Again,  rotating  the  pedestal  by  104.4  degrees  deepens 
the  null  at  130  degrees  to  about  48  dB  while  conserving  the  remaining  reductions  of 
about  11  dB,  20  dB,  and  11  dB  for  the  50,  90,  and  170  degree  nulls.  As  previously 
noted,  the  bistatic  scattering  widths  presented  Figure  11  demonstrate  how  a  tapered 
impedance  sheet  not  only  introduces  controllable  and  particular  nulls  in  the  bistatic 
scattering  widths,  but  the  ability  to  significantly  increase  the  reduction  amount  at 
certain  angles. 

4.1.2  PEC  pedestal  results. 

Consider  a  PEC  pedestal  having  a  radius  op  =  A  surface  gap-wrapped  (i.e.,  op  — 
op  =  .25A)  with  the  tapered  impedance  sheet  distribution  Zs((j))  shown  in  Figure 
9.  As  will  be  shown  below  for  this  more  obtrusive  pedestal,  the  relatively  simple, 
physically-realizable  complex  impedance  sheet  distribution  can  effectively  introduce 
specific  nulls  in  the  bistatic  scattering  width  unattainable  with  a  constant  impedance 
sheet  and  allows  the  null  positions  to  be  readily  controlled. 

The  bistatic  scattering  widths  illustrated  in  Figure  12  depict  the  TM2  polarization. 
At  its  initial  setup  state  the  single  sinusoid  impedance  sheet  distribution  introduces 
an  effective  null  near  55  degrees.  This  null  improves  the  two-dimensional  scattering 
width  by  about  16  dB.  As  before,  rotating  the  pedestal  by  36.0  degrees  deepens  the 
null  at  55  degrees  to  a  non-trivial  amount  of  40  dB.  Then,  rotating  the  pedestal  by 
108.0  degrees  effectively  repositions  the  null,  resulting  in  improvements  ranging  from 
about  80  to  180  degrees.  The  reduction  amount  increases  from  5  dB  at  80  degrees  to 
40  dB  at  180  degrees,  with  a  maximum  reduction  of  46  dB  at  170  degrees.  Again,  the 
bistatic  scattering  widths  presented  Figure  12  demonstrate  how  a  tapered  impedance 
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Figure  12.  2D  Bistatic  Scattering  Width  of  TM2  polarization  (Top:  0  degrees,  Middle: 
36.0  degrees,  Bottom:  108.0  degrees)  for  the  (^-dependent  impedance  sheet  partially 
filled  with  a  PEC  pedestal. 
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Figure  13.  2D  Bistatic  Scattering  Width  of  TEZ  polarization  (Top:  0  degrees,  Middle: 
90  degrees,  Bottom:  104.4  degrees)  for  the  (/(-dependent  impedance  sheet  partially  filled 
with  a  PEC  pedestal. 
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sheet  not  only  introduces  controllable  and  particular  nulls  in  the  bistatic  scattering 
widths,  but  the  ability  to  significantly  increase  the  reduction  amount  at  certain  angles. 

Figure  13  illustrates  the  bistatic  scattering  widths  of  TE^  polarization.  At  its  ini¬ 
tial  setup  state  the  single  sinusoid  impedance  sheet  distribution  introduces  a  general 
decrease  in  the  bistatic  scattering  width  by  at  least  5  dB  from  50  to  100  degrees,  with 
more  substantial  null  of  12  dB  centered  about  65  degrees  These  reductions,  however, 
come  as  a  trade-off  given  that  the  there  actually  exists  an  increase  in  the  bistatic 
scattering  width  from  20  to  40  degrees.  Rotating  the  pedestal  by  90  degrees  effec¬ 
tively  repositions  the  previous  null  and  introduces  additional  nulls  at  45,  80,  125, 
and  165  degrees,  with  corresponding  reductions  of  14  dB,  15  dB,  14  dB,  and  15  dB. 
Additionally,  at  this  rotation  a  general  improvement  in  the  bistatic  scattering  width 
exists  from  about  70  to  180  degrees.  Again,  rotating  the  pedestal  by  104.4  degrees 
deepens  the  null  at  165  degrees  to  about  42  dB  while  conserving  the  remaining  reduc¬ 
tions  of  about  18  dB,  13  dB,  and  12  dB  for  the  45,  80,  and  125  degree  nulls.  As  with 
every  figure  before,  the  bistatic  scattering  widths  presented  Figure  13  demonstrate 
how  a  tapered  impedance  sheet  not  only  introduces  controllable  and  particular  nulls 
in  the  bistatic  scattering  widths,  but  the  ability  to  significantly  increase  the  reduction 
amount  at  certain  angles. 

It  is  important  to  note  that  these  distributions  can  be  approximately  realized 
with  current  fabrication  techniques,  as  developed  in  [2],  The  employed  technique 
utilizes  a  constant  impedance  sheet,  typically  a  Kapton  KC  film,  and  manipulates  the 
effective-medium  properties  via  a  method  of  printing  metallic  squares  and  burning  out 
squares  of  various  dimensions  and  spatial  variations.  Though  the  presented  impedance 
sheet  distributions  are  physically  realizable,  the  fabrication  process  itself  has  some 
prohibitive  limitations  for  this  research.  First,  in  order  to  validate  these  2D  bistatic 
scattering  width  patterns,  a  2D-to-3D  conversion  must  be  employed  (in  the  absence 
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of  an  available  industry  computer  code),  as  shown  below 


2  l2 

03 D  ~  ~^—&2 Di 


(235) 


where  L  is  total  axial  length  of  the  cylinder  such  that  L  »  A.  Thus,  for  X-band 
frequencies  the  cylinder  must  have  a  minimum  length  of  36.5  cm,  as  determined  by  the 
minimum  frequency  of  8.2  GHz.  The  maximum  area  that  can  currently  be  fabricated, 
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Figure  14.  Example  of  a  tapered  impedance  sheet  based  on  a  Kapton  KC  film  substrate. 


however,  is  about  10  cm  sq.  Increasing  the  desired  frequency  range  to  circumvent  this 
limitation  is  not  viable  either  given  the  physical  limits  in  printing  microscopic  metallic 
squares  and  burning  out  microscopic  squares  in  the  Kapton  KC  him.  Despite  this 
current  limitation,  a  viable  electromagnetic  characterization  technique  can  still  be 
utilized  to  help  verify  the  currently  realizable  physical  samples. 


4.2  Electromagnetic  Characterization  Process 

In  order  to  mitigate  the  inhomogeneity  of  the  tapered  R-cards,  a  novel  electromag¬ 
netic  characterization  technique  must  be  adopted.  This  specialized  technique  employs 
a  stepped  flange  waveguide  and  modal  analysis  approach  in  order  to  approximately 
treat  the  tapered  R-card  as  a  piecewise  homogeneous  material.  The  main  objective  is 
to  focus  the  incident  electromagnetic  energy  into  a  smaller  sample  area  of  the  tapered 
R-cards  suspended  between  two  flanges,  as  shown  below. 
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Figure  15.  Theoretical  stepped  flange  waveguide  setup  where  the  gray  substance  rep¬ 
resents  the  tapered  R-card. 

This  process  ultimately  produces  the  theoretical  scattering  parameters  (Sf^ ,  S^y)  for 
the  entire  stepped  flange  waveguide  setup.  With  theoretical  values  for  the  complex 
permittivity  and  permeability  of  the  R-card  embedded  within  the  derived  scattering 
parameters,  an  iterative  root  search  must  then  be  implemented  to  solve  for  the  actual 
values  for  the  R-card’s  complex  permittivity  and  permeability.  This  is  accomplished 
via  the  Newton- Raphson  root  search  which  solves  for  the  following  relations 


e)  —  <S'iiP(a>)  <  tol 
-  Se2ip(u>)  <  tol 


(236) 


where  tol  =  10'.  With  the  modal  analysis  technique  the  fields  in  each  region  of  the 
stepped  waveguide  setup  are  expanded,  followed  by  satisfying  the  interfacial  boundary 
conditions  and  applying  an  appropriate  integral  testing  operator.  The  introduction  of 
the  stepped  flanges  disrupt  the  waveguide  wall  surface  currents  and  excite  higher-order 
modes.  These  higher-order  modes,  though  accommodated  by  the  mode-matching 
technique,  evanesce  within  the  waveguide.  This  can  lead  to  insufficient  transmitted 
energy  received  by  the  vector  network  analyzer  to  accurately  determine  the  complex 
permittivity  of  the  sample  R-card.  Two  solutions  to  this  problem  are  presented.  The 
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first  method  is  to  fill  the  stepped  flange  waveguide  with  a  low  loss  dielectric  ( e/j«,Ah) ) 
in  order  to  couple  more  energy  through  the  waveguide.  The  second  is  to  simply 
make  the  stepped  flanges  less  thick  (r),  such  that  the  energy  has  less  opportunity  to 
evanesce. 

The  ultimate  objective  is  to  determine  the  R-card  sheet  impedance.  This  can  be 
readily  accomplished  via  an  AC  expression  relating  the  effective  sheet  impedance  Zes 
and  the  extracted  complex  permittivity  eefj  [3] 


Zeff  — 


JVo  c 


aeSr  j u(eeS  -  e0)  ur(ef  -  1)  ’ 


(237) 


where  aefj  is  the  effective  complex  conductivity,  r  is  the  sheet  thickness,  u>  is  the 
angular  excitation  frequency,  c  is  the  speed  of  light,  and  =  eefj/eo  =  (e  +  a / jcn)/eo 
is  the  relative  effective  complex  permittivity. 


4.2.1  Experimental  Setup. 

All  measurements  were  conducted  with  a  HP8510C  Vector  Network  Analyzer  con¬ 
nected  to  a  WR90  waveguide  setup,  where  the  WR  represents  “waveguide  rectangu¬ 
lar”  and  90  indicates  the  aperture  width  of  0.90  inches.  The  full  X-band  frequency 
range  of  8.2  GHz  to  12.4  GHz  was  used  along  with  employing  the  industry  standard 
TRL  calibration  technique.  Furthermore,  all  measurements  are  accompanied  by  a 
differential  error  analysis  focusing  on  the  uncertainty  of  the  plate  flange  dimensions 
and  R.-card  thickness [22], 

Without  being  able  to  obtain  a  sample  tapered  R.-card,  a  proof  of  concept  was 
conducted  with  readily  available  uniform  R-cards  of  various  values.  Validation  of 
the  electromagnetic  characterization  technique  required  the  use  of  several  variations 
amongst  the  aperture  size  of  the  stepped  flange,  the  thickness  of  the  flange,  and  the  R.- 
card  resistivity.  Consequently,  each  particular  test  setup  variant  conducted  tests  with 
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Figure  16.  Experimental  setup  of  the  stepped  waveguide  and  flanges. 

three  different  R-cards  of  about  400  Q/sq,  900  Q/sq,  and  2000  fl/sq  to  demonstrate 
usefulness  across  a  range  of  impedances.  The  setup  variants  themselves  consisted  of 
plate  flanges  (6.47  mm  thick)  and  sheet  flanges  (.673  mm).  Furthermore,  both  the 
plate  and  sheet  flanges  had  three  more  aperture  size  variants  of  0.1  inches,  0.2  inches, 
and  0.3  inches.  Each  variant,  in  turn,  was  ultimately  compared  to  the  results  from  a 
full  X-band  flange  (t  =  0.9  inches).  With  the  full  X-band  flange  only  the  dominant 
TE10  mode  would  be  excited,  acting  as  an  ideal  baseline  measurement. 

4.2.2  0.1  inch  Aperture  Flange. 

The  0.1  inch  aperture  flange  represented  the  smallest  taper  to  be  experimented 
upon  and  performed  the  worse  of  the  three  aperture  widths.  With  such  a  small 
aperture  very  little  energy  is  transmitted  through  the  sample,  with  nearly  all  the 
energy  being  reflected  back,  as  illustrated  in  Figure  18.  At  this  width,  both  the  empty 
and  filled  plate  flanges  failed  to  allow  sufficient  energy  to  pass  through,  resulting  in 
the  S21  parameter  falling  completely  into  the  noise  floor.  The  sheet  flange,  however, 
allowed  sufficient  energy  to  be  transmitted  through  the  sample  to  measure  above  the 
noise  floor. 


67 


Figure  17.  Experimental  setup  variant  of  the  0.1  inch  plates  and  sheets. 
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Figure  18.  Experimental  scattering  parameters  for  the  0.1  inch  aperture  flange  (top- 
left:  empty  plate,  top-right:  filled  plate,  bottom-left:  sheet,  bottom-right:  X-band). 


Due  to  the  excessive  noise  present  in  the  scattering  parameters  from  the  0.1  inch 
aperture  flange  plates,  both  filled  and  empty,  it  is  impossible  to  employ  the  numerical 


Newton- Raphson  root  search  technique.  The  root  search  can,  however,  be  applied 
to  the  0.1  inch  sheet  flange  scattering  parameters,  to  extract  the  complex  impedance 
values  in  an  erratic  fashion.  These  extracted  impedance  values  fail  to  converge  to  the 
baseline  X-band  impedance  values,  as  evidenced  by  Figures  19,  20,  and  21. 


Figure  19.  Real  and  imaginary  components  of  the  effective  impedance  for  the  400  ohm 
R-card  with  a  0.1  inch  aperture  flange. 
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Figure  20.  Real  and  imaginary  components  of  the  effective  impedance  for  the  900  ohm 
R-card  with  a  0.1  inch  aperture  flange. 
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Figure  21.  Real  and  imaginary  components  of  the  effective  impedance  for  the  2000 
ohm  R-card  with  a  0.1  inch  aperture  flange. 

A  summary  of  the  maximum  difference  between  the  0.1  inch  sheet  flange  and  full 
X-band  extracted  complex  impedance  values  is  shown  below. 

Table  1.  Summary  of  Maximum  Complex  Impedance  Error  from  the  0.1  inch  Sheet 
Flange 


R-card 

Real  Component 

Imaginary  Component 

400  n 

27.5% 

11.4% 

900  n 

57.0% 

25.0% 

2000  n 

137.5% 

52.7% 

For  the  sheet  flange  measurements,  the  relatively  minor  noise  present  in  the  S21 
parameters  with  magnitudes  on  the  order  of  10~4  prohibit  the  root  search  technique 
from  effectively  extracting  the  correct  complex  impedance  values. 

4.2.3  0.2  inch  Aperture  Flange. 

The  0.2  inch  aperture  flange  represented  the  second  taper  to  be  experimented 
upon  and  performed  with  mixed  success.  For  all  flange  variants  sufficient  energy 
passed  through  the  sample  to  measure  above  the  noise  floor,  as  illustrated  in  Figure 
23.  In  this  case,  the  experimentally  measured  S21  magnitudes  for  the  aperture  flange 
plates,  both  filled  and  empty,  contain  relatively  significant  fluctuations  due  to  noise. 
The  sheet  flange  S21  magnitudes,  though,  are  seen  to  be  free  of  any  such  fluctuations. 
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Figure  22.  Experimental  setup  variant  of  the  0.2  inch  plates  and  sheets. 


Figure  23.  Experimental  scattering  parameters  for  the  0.2  inch  aperture  flange  (top- 
left:  empty  plate,  top-right:  filled  plate,  bottom-left:  sheet,  bottom-right:  X-band). 
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Also,  close  examination  upon  the  S21  magnitude  of  the  empty  to  filled  plate  flange 
reveals  that  the  included  plexiglass  allows  about  10  times  more  energy  to  transmit 
through  the  sample. 

Due  to  the  noise  present  in  the  scattering  parameters  from  the  0.2  inch  aperture 
flange  plates,  both  filled  and  empty,  the  root  search  technique  performs  unsuccessfully, 
as  previously  experienced  by  the  0.1  inch  sheet  flange.  Still,  the  dielectric  filled  plate 
flange  shows  an  improvement  in  the  extracted  complex  permittivity  over  the  empty 
plate  flange,  due  to  the  additional  energy  transmitted  through  the  sample.  The 
0.2  inch  sheet  flange  S21  parameters,  absent  of  the  noise  fluctuations,  effectively 
produces  extracted  complex  permittivity  values  that  converge  to  the  baseline  X-band 
permittivity  values,  as  evidenced  by  Figures  24,  25,  and  26. 


Figure  24.  Real  and  imaginary  components  of  the  effective  impedance  for  the  400  ohm 
R-card  with  a  0.2  inch  aperture  flange. 
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Figure  25.  Real  and  imaginary  components  of  the  effective  impedance  for  the  900  ohm 
R-card  with  a  0.2  inch  aperture  flange. 


Figure  26.  Real  and  imaginary  components  of  the  effective  impedance  for  the  2000 
ohm  R-card  with  a  0.2  inch  aperture  flange. 


A  summary  of  the  maximum  difference  between  the  0.2  inch  empty  plate,  filled 
plate,  and  sheet  flange  to  the  full  X-band  extracted  complex  impedance  values  is 
shown  below. 

Table  2.  Summary  of  Maximum  Complex  Impedance  Error  from  the  0.2  inch  Empty 
Plate  Flange 


R-card 

Real  Component 

Imaginary  Component 

400  n 

77.6% 

248.2% 

900  n 

205.7% 

403.2% 

2000  n 

553.8% 

890.5% 

For  the  empty  and  filled  plate  flange  S21  parameters,  the  noise  present  in  the  mea¬ 
surements  prohibit  the  root  search  technique  from  effectively  extracting  the  correct 
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Table  3.  Summary  of  Maximum  Complex  Impedance  Error  from  the  0.2  inch  Filled 
Plate  Flange 


R-card 

Real  Component 

Imaginary  Component 

400  ft 

51.9% 

173.0% 

900  ft 

128.4% 

293.5% 

2000  ft 

403.4% 

617.4% 

Table  4.  Summary  of  Maximum  Complex  Impedance  Error  from  the  0.2  inch  Sheet 
Flange 


R-card 

Real  Component 

Imaginary  Component 

400  ft 

6.5% 

11.7% 

900  ft 

11.9% 

12.1% 

2000  ft 

27.5% 

23.8% 

complex  impedance  values.  The  sheet  flange  S21  parameters,  without  any  significant 
noise,  produce  fairly  accurate  real  impedance  values  for  each  R-card  sample.  The 
imaginary  component,  on  the  other  hand,  tends  to  converge  to  the  incorrect  value  as 
it  is  significantly  more  sensitive  to  slight  fluctuations  in  the  measurements. 


4.2.4  0.3  inch  Aperture  Flange. 

The  0.3  inch  aperture  flange  represented  the  largest  taper  to  be  experimented 
upon  and  performed  with  the  most  success. 


Figure  27.  Experimental  setup  variant  of  the  0.3  inch  plates  and  sheets. 
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For  all  flange  variants,  with  the  exception  of  the  empty  plate,  sufficient  energy  passed 
through  the  sample  to  measure  above  the  noise  floor,  as  illustrated  in  Figure  28.  In 
this  case,  none  of  the  experimentally  measured  S21  magnitudes  appear  to  contain  sig¬ 
nificant  fluctuations  due  to  noise.  As  previously  seen,  comparing  the  S21  magnitude 
of  the  empty  to  filled  plate  flange  reveals  that  the  included  plexiglass  allows  about 
10  times  more  energy  to  transmit  through  the  sample. 


Figure  28.  Experimental  scattering  parameters  for  the  0.3  inch  aperture  flange  (top- 
left:  empty  plate,  top-right:  filled  plate,  bottom-left:  sheet,  bottom-right:  X-band). 


Without  any  apparent  noise  in  the  S21  parameters  from  the  filled  plate  and  sheet 
variants,  it  can  be  expected  that  root  search  technique  would  perform  successfully  for 
all  cases.  The  extracted  complex  impedance  from  the  empty  plate  flange,  however, 
fails  to  converge  to  the  correct  values,  as  evidenced  by  Figures  29,  30,  and  31.  This 
is  again  due  to  insufficient  energy  coupling  through  the  sample. 
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Figure  29.  Real  and  imaginary  components  of  the  effective  impedance  for  the  400  ohm 
R-card  with  a  0.3  inch  aperture  flange. 


Figure  30.  Real  and  imaginary  components  of  the  effective  impedance  for  the  900  ohm 
R-card  with  a  0.3  inch  aperture  flange. 


Figure  31.  Real  and  imaginary  components  of  the  effective  impedance  for  the  2000 
ohm  R-card  with  a  0.3  inch  aperture  flange. 
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A  summary  of  the  maximum  difference  between  the  0.3  inch  empty  plate,  filled 
plate,  and  sheet  flange  to  the  full  X-band  extracted  complex  impedance  values  is 
shown  below. 

Table  5.  Summary  of  Maximum  Complex  Impedance  Error  from  the  0.3  inch  Empty 
Plate  Flange 


R-card 

Real  Component 

Imaginary  Component 

400  n 

15.4% 

75.2% 

900  n 

51.2% 

124.6% 

2000  0 

167.9% 

280.0% 

Table  6.  Summary  of  Maximum  Complex  Impedance  Error  from  the  0.3  inch  Filled 
Plate  Flange 


R-card 

Real  Component 

Imaginary  Component 

400  n 

5.2% 

4.6% 

900  n 

3.2% 

4.3% 

2000  n 

8.7% 

8.1% 

Table  7.  Summary  of  Maximum  Complex  Impedance  Error  from  the  0.3  inch  Sheet 
Flange 


R-card 

Real  Component 

Imaginary  Component 

400  fl 

3.8% 

1.9% 

900  n 

6.8% 

3.6% 

2000  n 

19.6% 

9.2% 

For  the  empty  plate  flange  S21  parameters,  the  inadequate  amount  of  transmitted 
energy  prohibits  the  root  search  technique  from  effectively  extracting  the  correct 
complex  permittivity  values.  The  hlled  plate  and  sheet  flange  S21  parameters,  with 
sufficient  energy  transmitted,  produce  very  accurate  complex  permittivity  values  for 
each  R-card  sample.  In  each  case  where  the  complex  impedance  extraction  algorithm 
failed,  the  experimentally  measured  S21  values  fell  near  0.1  or  below.  Thus,  the  exper¬ 
imental  results  collectively  confirm  a  material  measurement  principle  which  requires 
that  at  least  10%  of  the  incident  electromagnetic  energy  be  transmitted  through  the 
test  sample  for  precision  measurements. 
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4.2.5  Summary. 


While  the  Newton- Raphson  root  search  is  extremely  accurate,  it  requires  an  ini¬ 
tial  guess  of  the  unknown  parameter.  Without  prior  knowledge  of  the  material  under 
test  the  technique  may  not  converge  or  could  even  converge  to  the  incorrect  solution. 
This  problem  is  exacerbated  when  using  relatively  small  scattering  parameter  mea¬ 
surements.  With  relatively  small  values  the  amount  of  uncertainty  present  has  a  more 
substantial  effect  on  the  root  search  algorithm,  causing  it  to  become  unstable.  Addi¬ 
tionally,  a  minimum  S21  measurement  of  0.1  is  required  for  precision  electromagnetic 
characterization.  Filling  the  stepped  flange  waveguide  with  a  low  loss  dielectric,  in 
this  case  plexiglass  (er  =  2.6),  allows  about  ten  times  as  much  incident  energy  to 
couple  through  the  sample.  Simply  decreasing  the  thickness  of  the  stepped  flange, 
however,  is  more  effective  in  permitting  energy  to  couple  through.  The  difference  can 
be  expected  by  considering  the  forward  propagation  component  of  the  electromagnetic 
wave, 


e  lmZ  for  m  —  1, . . . ,  M  (238) 

r)m  \J ^‘xrn  ^  ^  ~  \J j  kXm  (2771  l)~j  (239) 

where  c  is  the  speed  of  light,  er  is  the  relative  permittivity  of  plexiglass,  fir  is  the 
relative  permeability  of  plexiglass,  and  t  is  the  thickness  the  waveguide  width.  Com¬ 
puting  the  amount  of  attenuation  for  the  first  m  =  1, . . . ,  10  modes  when  passed 
through  the  plate  and  sheet  flanges  produces  the  following  figures. 
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Figure  32.  Percentage  of  first  10  modes  to  pass  through  the  0.3  inch  stepped  flange 
(top-left:  empty  plate  flange,  top-right:  filled  plate  flange,  bottom:  sheet  flange. 


Analyzing  Figure  32  confirms  that  simply  decreasing  the  flange  thickness  more  signif¬ 
icantly  increases  the  amount  of  electromagnetic  energy  coupled  through  the  sample. 
For  both  the  empty  and  filled  flange  cases,  only  the  first  two  modes  succeed  in  cou¬ 
pling  through  the  sample,  with  the  the  second  mode  transmitting  only  at  0.04.  Yet, 
the  sheet  flange  effectively  allow  an  appreciable  amount  of  energy  to  be  transmitted 
through  nearly  all  ten  modes. 
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V.  Conclusions  and  Recommendations 


5.1  Conclusions 

This  research  investigated  the  ability  to  simultaneously  alter  the  real  and  imag¬ 
inary  components  of  an  impedance  sheet  to  achieve  novel  distributions  for  precision 
scattering  and  shielding  control  and  manipulation.  The  theoretical  development  de¬ 
pended  on  an  analytic  approach  for  calculating  the  scattering  characteristics  from 
cylindrical  sheet  impedances  having  a  0  dependency  with  free-space  filled,  partially- 
filled  PEC  filled,  and  general  dielectric  filled  cases.  It  relied  on  enforcement  of  general 
impedance  sheet  boundary  conditions  in  addition  to  a  Fourier  series  representation 
of  the  impedance  sheet  distribution  with  Wronskian  and  recurrence  relations.  For 
electromagnetic  characterization,  the  additional  degree  of  complexity  caused  by  the 
(f)  dependency  was  successfully  mitigated  in  a  piecewise  homogeneous  approach  via 
a  reduced  aperture  waveguide  system.  This  stepped  waveguide  technique  utilized  a 
modal  analysis  to  account  for  the  excited  higher-order  modes,  effectively  allowing 
smaller  cross-sections  of  the  ^-varying  impedance  sheet  to  be  characterized.  Further¬ 
more,  it  was  established  that  10%  of  the  incident  electromagnetic  energy  at  minimum 
must  transmit  through  the  sample  under  test  to  achieve  precision  electromagnetic 
characterization. 

5.2  Future  Research 

For  the  novel  impedance  sheet  distributions,  advances  in  fabrication  techniques 
will  eventually  allow  larger  tapered  R-card  samples  to  be  made.  Then,  the  scattering 
width  2D-to-3D  conversion  may  be  employed  given  that  the  minimum  length  crite¬ 
rion  of  35  cm  for  X-band  frequencies  has  been  met  and  experimental  testing  may  be 
completed.  Another  point  of  interest  for  the  impedance  sheet  distributions  is  the  idea 
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of  multiple  layers  to  possibly  achieve  less  frequency  dependence.  This  setup  would 
required  a  layered  Green’s  function  analysis  with  the  resistive  sheet.  Furthermore, 
it  may  be  possible  to  achieve  an  analytic  expression  to  relate  both  the  position  and 
reduction  of  effective  nulls  in  the  two-dimensional  bistatic  scattering  width  if  the 
general  approach  is  limited  in  scope  to  a  single  sinusoid.  With  the  material  charac¬ 
terization  aspect,  the  presented  reduced  aperture  waveguide  system  could  possibly  be 
improved  upon  by  allowing  the  Vector  Network  Analyzer  to  average  the  experimental 
scattering  parameters  over  100  or  more  measurements.  An  additional  material  char¬ 
acterization  technique  that  could  be  explored  involves  the  tapered  impedance  sheet 
being  rotated  180  degrees  in  a  typical  EMC/EMI  setup  while  monostatic  measure¬ 
ments  are  recorded  for  every  increment  of  rotation.  This  method  would  require  an 
adaptation  of  a  method  of  moments  approach  to  characterize  the  discertized  tapered 
impedance  sheet. 
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Appendix  A.  Computationally-Efficient  Formulation 


The  cylindrical  wave  transformation  representation  of  the  electromagnetic  fields, 
TM2  and  TE2,  can  be  written  as 


E‘(p)  =  zEle-i^  =  zEl  ]T 


(240) 


H‘(p)  =  =  zHl  Y, 


(241) 


which  can  be  manipulated  into  a  more  computationally  efficient  form  through  the 
use  of  certain  Bessel  function  properties.  Expanding  the  cylindrical  wave  function 
representation  into  three  terms  of  negative  indices,  zero  index,  and  positive  indices 
results  in 


E  rnJn(Mejn0  =  Jo(M  +  E  rnUk0p)ejnt  (242) 


Re-indexing  the  last  term  for  n  — >  —n  produces 


E  rnUkop)ejn *  =  EfJ-»(M< 


(243) 


Ejn(-l)nJn(Me'^ 


(244) 


(245) 


hn  »(M< 


(246) 
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Substituting  (246)  back  into  (242) 


5]  rnUk0p)ejn *  =  J0(M  +  Y  rnUkop)(ejn<t>  +  e 

7.— — oo  n=  1 

oo 

=  Jo(M  +  Y.  j~nJn(fcop)(2cosn0) 

n—  1 
oo 

=  Y  enj~nJn(k0p)  COS  710, 


~jn4 


n= 0 


where 


(247) 

(248) 

(249) 


{1  . . .  n  =  0 

2  ...n/  0 


(250) 


The  incident  electromagnetic  fields  can  now  be  represented  in  a  more  computationally 
efficient  form  as 


E\p)  =  zE'0ze-ikx  =  zEqz  Y  £nT  "  Jn(M  cos  n0,  (251) 

n= 0 
oo 

H\p)  =  zH'0ze~jkx  =  zH^Y^^nMcosruj).  (252) 

n= 0 
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Appendix  B.  Wronskian  and  Recurrence  Relations 

Developed  for  the  TM~  field  set,  the  following  will  be  additionally  simplified  by 
utilizing  Wronskian  and  recurrence  relations,  such  that 

J n(k0a)RW(k0a)  -  J'n(k0a)R^ (k0a) 

=  J n(k0a)  [J 'n(k0a)  -  jY'n(k0a )]  -  J 'n(k0a)  [ Jn(k0a )  -  jYn(k0a)\ ,  (253) 

=  Jn(k0a)J'n(k0a)  -  jJn(k0a)Y'n(k0a)  -  yn(k0a)Jn(k0a)  +  jJ'n(k0a)Yn(k0a),  (254) 
=  j  [J'n(k0a)Yn(k0a)  -  Jn(k0a)Y'n(k0a )] ,  (255) 

where  Hankel  functions  of  the  second  kind  and  their  derivatives  can  be  defined  as 

H £\k0a)  =  J n(k0a)  -jYn(k0a),  (256) 

H i2)'(M  =  yn(k0a)  - ]Y'n(k0a).  (257) 

Applying  the  following  derivative  Bessel  function  identity, 

Tl 

Bn(^)  =  -Bn+i(z)  +  -B  n(z),  (258) 
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to  (255)  results  in 


jn(MH^'(M  -  J  UMHi2)(M 


=  3 


n 


-Jn+i (k0a)  Y  — J n(k0a)  )  Y n(k0a) 


n 


-J n(k0a)  (  -Yn+i (k0a)  Y  —  Y n(k0a) 


=  3 


n 


-Jn+i(/coa)Yn(/c0a)  +  - —  Jn(k0a)Yn(k0a) 

Koa 


n 


Y Jn (fcoo) Yn_)_i  (k§(i)  Y  .  Jn (^-0®) Yn. (^0®) 

k^a 

=  3  [Jn(k0a)Yn+1(k0a)  -  Jn+i(k0a)Yn(k0a )] , 

=  J(-1)  2 


Tik^a 


(259) 


(260) 

(261) 

(262) 


Thus,  utilizing  both  Wronskian  and  recurrence  relations  leads  ultimately  to  the  fol¬ 
lowing  relation 


jn(MHi2)'(M  -  JUMHi2)(M  = 


(263) 


Developed  for  the  TE2  field  set,  the  following  will  be  additionally  simplified  by 
utilizing  Wronskian  and  recurrence  relations,  where 


4(MhL2)(M  -  J„(MHl2)'(M 

=  yn(k0a)  [J n(k0a)  -  jYn(k0a )]  -  Jn(&oa)  [J^&oa)  -  jY'J&oa)] ,  (264) 

=  J(,(/c0a)Jn(^oa)  -  jJn(^oa)Yn(^0a)  -  Jn(k0a)J'n(k0a)  Y  jJn(k0a)Y'n(k0a),  (265) 
=  3  [Jn(k0a)Y'n(k0a)  -  J'n(k0a)Yn(k0a)J ,  (266) 
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where  Hankel  functions  of  the  second  kind  and  their  derivatives  can  be  defined  as 


H £\k0a)  =  J n(k0a)  -jYn(k0a), 
Hi2)'(M  =  KM  -]Y'n{k0a). 


(267) 

(268) 


Applying  the  following  derivative  Bessel  function  identity, 


Tl 

Bn(^)  =  -Bn+i(z)  +  -B  n(z), 


(269) 


to  (266)  results  in 


yn(k0a)R^\k0a)  -  Jn(k0a)H.W(k0a) 
—  J 


n 


J  n(k0a)  -Y  n+i(k0a)  +  — Y  n(k0a) 

koa 


n 


"Jn+l(^0®)  "h  ,  Jn(^0®)  )  Yn(koa) 

koa 


=  J 


n 


■  Jn(/CoCl) Yn_)_i  (/uoQ)  Y  .  Jn (A/0®) An (^0®) 

koa 


n 


Y  Jn+l  (A/0®)  An(^o®)  ,  Jn(^0®)An(^0®) 

koa 

J  [  Jn+l  ( koa)  Yn  (/coq)  Jr)  (^0®)  Yn+1  (/coo)]  , 

2 

=  j 


irkoa 


(270) 


(271) 

(272) 

(273) 


Thus,  utilizing  both  Wronskian  and  recurrence  relations  leads  ultimately  to  the  fol¬ 
lowing  relation 


j'n(MHi2)(M  -  Jn(M42),(M  = 


(274) 
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